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Abstract 

We develop the phenomenological amplitude of the vrA^ ttttN reaction describing the 
exchanges of A and N^, along with the OPE mechanism. The contribution of the latter con- 
tains 4 independent low energy parameters (up to 0(/c*^) order). The terms of the polynomial 
background are added to stand for far resonances and for contact terms originating from the 
off-mass-shell interactions. These terms are introduced with the account of isotopic, cross- 
ing, C, P and T symmetries of strong interactions. 

The data consisting of total cross sections in the energy region 0.300 < PLab ^ 500 MeV/c 
and ID distributions from the bubble-chamber experiments for three reaction channels were 
undergoing fittings to determine free parameters of the amplitude. The best solutions are 
characterized by Xbf ~ 1-16. At the considered energies the isobar exchanges are found to 
be more important than OPE. The obtained solutions reveal the need in more precise data 
and/or in polarization measurements because of large correlations of isobar parameters with 
the OPE ones. 

The theoretical solutions were used for modeling the Chew-Low extrapolation and the 
Olsson-Turner threshold approach. It is shown that the noncritical application of the former 
results in 100% theoretical errors, the extracted values being in fact the random numbers. 
The results of the Olsson-Turner method are characterized by significant systematic errors 
coming from unknown details of isobar physics. 
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1 Introduction 



The pion-production reactions play the important role in the low energy physics of ele- 
mentary particles and nuclei. The elementary processes vriV — > imN, 'yN — > vrvrA^ of pion 
production on nucleons are currently in the focus of investigations due to the progress of 
Chiral perturbation theory (ChPT) the foundation of which was created by the Weinberg's 



work ||81| and principal steps were done in the series of papers by Gasser and Leutwyler . 
(For more references and reviews of this approach as well as more recent trends one can use 
the book by Donoghue, Golowich an Holstein p9[ and the review papers by Meissner 



and Pich a deep insight into the contemporary interpretation of ChPT and the role of 
spontaneous symmetry breaking is provided by a series of preprints by Leutwyler [0.) 
The predictions of ChPT |^ (summary of the most interesting predictions as well as 



forthcoming experimental tests might be found in the talk |]68[) explain the very sharp interest 



in the values of vrvr-scattering lengths and other parameters of the vrvr scattering. Since it 
is not possible experimentally to create the pionic target or the colliding pion beams there 
are only indirect ways for obtaining experimental data on the vrvr scattering. The reactions 
hN —>■ TTTiN and K ttttcz/ are considered as the most important sources of the (indirect) 
information on low energy characteristics of the vrvr interaction. (The former reaction in 
what follows will be simply referred as 7r27r when possible.) 

The review of methods of extracting the latter characteristics from 7r27r data being in 
use previously and details of their applications might be found in the paper |50| by Leksin. 
The methods considered for application to the modern experiments are the following . 



1. The Chew-Low extrapolation procedure by Goebel, Chew and Low is an appar- 
ently model-independent approach. It can provide the complete information on the vrvr cross 
section provided the OPE dominates and the interval of the nucleon momentum transfer r 
(which equals the mass of the virtual pion) allows an unique extrapolation. The last condi- 
tion and the need for sufficient statistics shifts the region of application of the Chew-Low 
procedure to rather large values of energy {^1-4 GeV). When comparing the phase space 
of momentum transfer — 20/i^ < r < — 0.2/i^ at pLab = 500 MeV/c with the distance of 
extrapolation ^ /i^ it becomes obvious that provided there are enough statistics, the kine- 
matics of the vr2vr reaction itself does not prevent the use of the Chew-Low procedure at 
much smaller energies. It is the presence of contributions like that of A and A^* isobars which 
makes a straightforward extrapolation difficult at moderate energies due to the perturbation 
of the simple r-dependence of the OPE graph. The absolute values of all other contributions 
are killed at the extrapolation point r = /i^, but the result of extrapolation is known to be 
sensitive to the shape of the extrapolating curve [0, |5^ . 

2. In view of importance of concurrent mechanisms at intermediate energies the approach 
might be changed to determining the OPE parameters directly in the physical region of the 
reaction — this was implemented by the model of Oset and Vicente- Vacas It is clear 
that the neglect of a specific resonance contribution and/or the account of another one are 
capable to provide a lot of derivatives of the Oset- Vicente model. 

There is the energy region below PLab = 500 MeV/c, where the variation of r is sufficient 
to detect the OPE contribution since the contributions of the concurrent processes (being 
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nonnegligent) are smooth enough. The model |]2T| takes these features into account and 
naturally completes the Oset- Vicente approach in this specific energy domain. 

3. The investigations by Olsson and Turner are confined to the threshold of 7r27r 
reactions. Since the phase space of r variable shrinks to the point tq = — 2.31/i^ the applica- 
tion of the Chew-Low procedure is impossible there. The idea is to take advantage of Chiral 
Dynamics at the tt2tt threshold. 

The important results of the approach are the formulae expressing the vrvr-scattering 
lengths in terms of the threshold characteristics of the pion-production reactions. These 
formulae have gained a broad scale of application, especially in the recent years when new 
data on the nN irirN reactions in the close-to-threshold energy region became available 

m. 



41, 43, 44 , 72, 53 



The evidence of the importance of next-to-leading order terms of Chiral Lagrangian for 
the nN interaction and, in particular, for the 7r27r amplitude makes it necessary to 
modify the Olsson-Turner method. Recently the approach of heavy baryon approximation 
was used to derive corrections to the Olsson-Turner formulae and to make direct predictions 
of ChPT for the threshold 7r27r amplitude itself |]TT |. 



The recent results [^] of the so called Generalized ChPT approach and the progress 



in the two-loop ChPT calculations [|14|, ^ are claiming for more precise experimental infor- 
mation on the TTTT interaction at low energies at 0{k^) order. Some experiments listed in ref. 
pE| had already been finished: BNL and LAMPF results on total cross sections have been 
published [Q, ||6^, ID-distributions have appeared recently in the WWW (home pages 
http:/ /helena. phys. Virginia. edu/~pipin/El 179/E1179. htm\ 



http:/ /helena.phys. Virginia, edu/ ~pipin/E857/E857. html ) , 



higher distributions are in progress, the off-line treatment of experimental tapes of the TRI- 
UMF experiment |]75| will be completed soon. Therefore, it is timely to make the solution 
which theoretical method from the above list is capable to provide more reliable treatment 
of the modern 7r27r experiments. 

The main goal of the present work is to provide grounds for comparison of the listed 
approaches. For this purpose we develop the most extensive phenomenological amplitude of 
the considered reaction suitable for near-threshold and intermediate energy regions. 

We are basing on the approach 2. since both the Chew-Low extrapolation and the 
Olsson-Turner threshold formulae can not provide any hint for cross-checking of the rest 
methods. We try to fix the phenomenological amplitude by fitting the data on total cross 
sections and distributions of the reaction in question in the energy region from threshold up 
to PLab < 500 MeV/c. 

To avoid any doubt in the results in respect to correctness of acceptances, systematic 
errors, etc. the distribution data are chosen to be the bubble-chamber ones. This leaves us 
with rather old experiments (which are discussed in sect. 3). However, the significant part 
of the data has never been published and most its part eluded strong theoretical analysis 
(apart authors' checks of some isobar-like models). Therefore, it seems important to develop 
the tools of theoretical treatment of such data for determination of characteristics of pion- 
pion and pion-nucleon interactions along with other parameters of the phenomenological 
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amplitude. 

It is worth noting that the considered reaction via unitarity relations is directly connected 
to the fundamental process of elastic pion-nucleon scattering at intermediate energies and 
other processes like 'fN — > ttttN which gain the raising interest in the ChPT approach. It 
also enters the description of pion-nuclei scattering as an elementary process. Therefore, the 
structure of 7r27r amplitude is of great importance for nuclear and particle physics. 

The paper is organized as follows. The content of sect. 2 reminds the basics of the 
low energy phenomenology of the vrA^ ttuN reaction and describes the structure of our 
amplitude. Sect. 3 provides the summary of experimental data on distributions and total 
cross sections which we are analyzing. Sect. 4 is devoted to specifics of the fitting procedure 
and main results of analysis. The results of modeling the Olsson-Turner and Chew-Low 
approaches are exposed in sects. 5, 6 along with the discussion of some properties of our am- 
plitude and theoretical solutions. The summary, the concluding remarks and the discussion 
of the perspectives of the further development are given in Conclusions. 



4 



2 Model of TiN -> titiN Amplitude 



The principal features of the near-threshold phenomenology of the nN -kttN reaction 
had been already discussed in the paper in the quoted work the smooth background 



amplitude (+ the OPE one) had been derived for the energy domain bounded by the reaction 
threshold and the threshold PLab ~ 500 MeV/c of the A-isobar creation in the final state. 
However, the statistically significant data on 7r27r distributions (described in the next section) 
exist just for the boundary (and slightly above) of the pointed energy region. This makes 
necessary to modify the amplitude elaborated in ref. since the smoothness assumption 
is hardly to be valid there. 

In the current section we recall the phenomenology of 712-7? processes to be taken into 
account (subsect. 2.2) and principal parts of the modified amplitude (subsect. 2.3) the 
parameters of which must be determined from the data fittings. We start with the brief 
description of the spin-isospin structure of the discussed amplitude. 



2.1 7tN ttttN Amplitude 
2.1.1 General Structure 

In the isotopic space the amplitude M^^^Xf] Aj) of the reaction 

7i%h) + N^{p;\,) ^ 7r^(A;2) + 7r^(A;3) + iV^(g;A/) (1) 

has 4 degrees of freedom and might be expressed either in terms of the definite isospin 
amplitudes or in terms of the isoscalar ones. Separating the nucleon spinor wave functions 



(2) 



where the (^75) multiplier ensures the correct P-parity properties of the considered amplitude 
one can define the isoscalar amplitudes A, B, C, D by 



m; 



abc 

13a 



AtIJ'^ + Pr^,5"^ + CtIX + D«e'^''=(5, 



/3a 



(3) 



The analysis |jT9| of the spinor properties of the amplitude (H) allows to express each of 
the isoscalar functions A, P, C, D in terms of 4 independent form factors in the crossing- 
covariant way 

A = SA + VAk + VAk + i/2TA[k,k] 
B = SB + VBk + VBk + i/2TB[k,k] 
C = Sc + Vck + Vck + i/2Tc[k,k] 
t) = SD + VDk + VDk + i/2TD[k,k] ■ (4) 
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Here, k, k are the crossing-covariant combinations of pion momenta 

k = -ki + ek2 + eks; k = -ki + ek2 + ek^ , (5) 

where e = exp(27ri/3) = —1/2 + i\/3/2, e = e* = —1/2 — i\/3/2. These combinations 
together with the independent crossing-invariant ones 

Q = -ki + k2 + k3=p-q; P=p + q (6) 

are used to define 5 independent crossing-covariant scalar variables 

jy = p.k;u = p.k, (7) 

which completely determine the point in the phase space of the considered reaction. The 
expressions of all scalar products of particles' momenta are given in the paper |^ for the case 
of the unbroken isotopic symmetry. The definitions (|^) , (P) , are assumed for the physical 
particles (^2 is the vr" momentum in the reactions { — \-n} and {—Op}', in the {+0p} case 
is the Tr"*" momentum). All actual kinematical calculations in the computer programs are 
being processed with the isotopic symmetry breaking due to the particles' masses — this 
complicates the expressions given in the quoted paper. For simplicity we shall hold on the 
unbroken isotopic symmetry case in the illustrations and in the discussions which follow. 

The r variable coincides with the mass of the virtual pion of the OPE graph. The Air 
vertex of this graph is characterized also by the Mandelstam variables. The discussion of 
the off-shell dependence of the Air vertex on these variables is given in the paper |^. To 
avoid ambiguity we use only the dipion invariant mass 



{k2 + h 



\2 



in the discussion below. 

The amplitudes of the observable channels of the reactions vrA^ —>■ ttttN with the conven- 
tion for the normalization of the particle states adopted in |21| are provided by relations (the 
nontrivial statistical factors are taken into account in the cases of the {OOn} and {+ + n} 
channels) : 



M{_+„| = V2/2 (A + C); M{oo„} = 1/2 (A) ; M{++„| = 1/2 (B + C) ; 

M{_op} = 1/2 (C - 2D) ; M^+op} = 1/2 {C + 2D) . (9) 

They have the same form as in eqs. (§) 

Mx = Sx+Vxk+Vxk+t/2Tx[k,k]; X = {{-+n}, {-Op}, {OOn}, {++n}, {+0p}) , (10) 

where the spinor structures are defined according to the expansions (|), e.g. = 
\^/2{Sa + Sc) , etc. In practice, the following combinations of the vector structures 

V^ = {Vx + Vx)/2; V^ = {Vx-Vx)/{2z) (11) 

are being used in the course of calculations. 
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2.1.2 Cross Section 



The experimental data of the channel X are compared with the theoretical cross sections 



0"^) for a given experimental point (a) 



d^ko 



d^q 



{27t)^5\p + ki-q-k2-ks) \\M\\'Q, 



AjJ (27r)32A;2o (27r)32A;3o (27r)32go 
Here, = {ficf = 0.38937966(23) [GeV^ mbarn] is the conversion constant, 



(12) 



4J = AJ{p-kiy 



stands for normalization of the initial state and the characteristic function B(q,) = B(q)(p, ki, 
Q, ^2, ^3) of the bin (a) describes the appropriate cuts in the phase space (if any; in the case of 
the total cross section this function is equal to 1). The notation | |M| p stands for the squared 
modulus of the amplitude summed over polarizations of the final nucleon and averaged over 
ones of the initial proton (we shall call it simply matrix element). The statistical factor 
(equal to the product of l/n^\ over subsets of identical particles) for calculation of the cross 
section was included into definitions (|^) of the physical reaction amplitudes. 

The matrix element ||M|p = ||Mx|P is the quadratic form of the vector of spinor struc- 
tures {Sx, Vx, Vx, Tx) (or, the same, of the vector {Sx, Vjf , Vj^, Tx))'- 

\\Mxf = l/2 5] [M(g;A/)Mx(^75)«(p;A0] [u(g; A/)Mx(^75)«(p; A.) 



(X = {- + n}, {-0 p}, {0 n}, {+ + n}, {+0 p}) ; (13) 
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(14) 



For the simplified case of equal pion masses the matrix G of the above form is given in 
the paper In practice, the matrix for every channel had been calculated separately with 
the physical masses of all particles in the considered channel. 

To plot the data and the theoretical results we define the quasi-amplitude (M(q,)) which 
is the square root of the cross section (|l^) divided by the phase space in the case of the 
total cross section data; in the case of distributions both the cross section and the phase 
space are independently normalized to 1 — we call this quantity normalized quasi-amplitude 

norm- 



(a) 



mi 



cr 



\ 



(a) 



(||M||^ 



(15) 



(a) I 



7 



Here, the phase space cr(l) is the theoretical cross section (|I2|) obtained with the unit matrix 
element. 



2.1.3 Threshold Properties 

At the threshold of the reaction there are considerable simplifications in the representations 
(D or (|l^) since: a) the momenta of the outgoing pions are equal to each other {k2 = k^), 
b) the contribution of the D amplitude to the amplitudes (^) becomes zero and c) B = C. 

Another simplification takes place in the sum over final polarizations and the average 
over initial ones of the amplitude (0) squared modulus. At the threshold it degenerates to 

\\Mxf = i-ro)[S'x - (2m + 3/i)(V^ + V^)]' , (16) 
(X = {- + n}, {-0 p}, {0 n}, {+ + n}, {+0 p}) , 

where Tq = — 3/i^m/(m + 2/i) is the threshold value of the variable r. This provides the 
grounds to introduce the threshold amplitudes. 
The isotopic threshold amplitudes are 

A° = ^[5° -(2m + 3^)(F° + V:4)] ; 

B' = V^oK - {2m + 3^^m + V^)] ■ (17) 
C = V^oK-i2m + 3fi){VS + VS)], 

where all form factors are to be calculated at the threshold values of kinematical variables 
— this provides 

5° = C° . (18) 



To link the isotopic threshold amplitudes {\17\) with the experimental information let us 
construct the quantities M° 

M°„+„j = {A^ + B')V2/2, Mf^ori = 572, Mfoo.} = ^72, 

M|7-,„} = {B' + C')/2 = B\ Ml,^^ = 572, (19) 

which include both the isotopic and the statistical factors. 

Then the absolute values of the threshold amplitudes might be expressed in terms 
of threshold limits of the experimental quasi-amplitudes (|T5|): 



(M|_+„})|_^ = \A' + i?°|/v^ , (M{-ori)|_„ = m/2 , (M{oon})|_„ = \A'\/2 , 

= |i?r (M{+ori)|_„ = |5°|/2. (20) 



The threshold amplitudes (p!9| ) are dimensional and in the following their numerical values 
will be given in [(GeV)"^]. 

The threshold limits (|T^, being determined by only two isotopic threshold amplitudes 
and 5°, must satisfy three relations. The first two are straightforward: 

\Ml^^}\ = 2|M°_o,}| = 2\Ml,^y\ . (21) 



8 



From the definition ( p!9D it follows that 



Mf^^^j = V2Ml_^^^ - 2M{'oo.} . (22) 

This implies that, depending on the relative sign of and 5°, the relation between positive 
quantities (^) might be either 

iMf^+^jl = v^|Mf_+„}| - 2|M{'oo„}| (23) 

or 

|M{V„j| = v^|MO_^„j| + 2|M{'oo„j| (24) 

or 

|Mf+_,„^| = -v^|M°_+„}| + 2|Mfoo„|| . (25) 

Now it is necessary to recall the known properties of 7r27r cross sections. Since in the near 
threshold region the known cross sections cTjoon} of the {0 n} channel are approximately 
equal to the cross sections cr{_+„} this excludes the possibility (^3j) because it results in 
the negative RHS. In the whole region pLab < 0.5GeV/c the cross sections o-j^^^n} are 
considerably smaller than cr{_+„,} and/or crjoon}; so, the variant ( p^ is impossible. Thus, 
combining with the previous two relations, we can state: 

2|M°_o^}| = 2|Mf_,o^^| = \Ml^^^\ = 2|M{Vj|-v^|M°_^„^| . (26) 



2.1 A Remarks on Sign Ambiguity 

To resume the discussion of the general structure of the 7r27r amplitude let us consider 
briefly the problem of the sign ambiguity in the theoretical amplitude which is used to fit the 
experimental cross section data. The examination of relations (EDf) shows that already in this 



simplified threshold case apart the overall sign ambiguity of isotopic threshold amplitudes 
y4°, B*^ their relative sign might also become indefinite depending on the accuracy of the 
experimental information. 

In what follows we call the solution physical (unphysical) when and are of different 
(equal) signs. In terms of the vrvr-scattering lengths the two cases of signs correspond to the 
different (equal) signs of ag^^ and Oq^^. 

In general there are no grounds to wait for a simplification to take place at a distance from 



the threshold. Indeed, the expression (|T^) for the matrix element ||Mx||^ might be brought 
to the diagonal form (for example, in terms of the analogues of the diagonal derivative 
amplitudes of Rebbi [^) in which it becomes the sum of four squared modules of the 
orthogonal amplitudes. The abundance of solutions found in the course of data fittings 
should be explained in part by a variety of choices of the signs in the above four terms in 
the cross section of every channel. 
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2.2 Near— Threshold irN irirN Phenomenology Guidelines 



1. The current knowledge of physics of the pion-nucleon interactions provides no evidence 
of possible mechanisms or processes resulting in the strong variation of the amplitude of 
the hN vrvrA^ reaction in the phase space at energies up to ^ 2m,r above threshold. All 
resonances but and A are located outside the phase space. 

2. The contribution of the near-threshold A pole is only due to the chain ttN — > A ^ 
nnN. It is suppressed 1) by the negligible width of the decay A nnN; 2) by Quantum 
mechanics: A is located precisely at the 7r27r threshold where the process tiN — > nnN must 
proceed through the waves Pu and P31 of the initial pion-nucleon system (in the L212J 
notation) while A can be created only in the P^-^ wave of the n-N system. 

3. The closest to the phase space contribution of the meson-resonance type is that of 
OPE — all the other are very distant at the discussed energies. 

4. The tails of resonances result in some constant background in the physical domain 
and, at most, in some slow variation of the amplitude. (One would need the extraordinary 
precision of experimental data to distinguish between a linear function in the corresponding 
variable and the far pole.) 

5. The number of free (or unknown) parameters in tttt, ttN, ttA, ttN^ interactions exceeds 
the number of degrees of freedom of the smooth near-threshold amplitude. 

6. All resonances contribute with terms built of various scalar products of momenta. 
There are at most five such products which might be considered to be independent or, the 
same, they are built of five invariant variables (at the fixed initial energy, four variables 
left). Hence, the arbitrary linear function is parametrized with 6 parameters: a constant + 
5 coefficients at linear variables. This applies to any of 4 independent isospin amplitudes 
which describe 5 observable channels of the reaction. 

7. There are four different spinor structures in the 7r27r amplitude. Even in the case of 
the unpolarized experiment they provide the specific dependence on variables which might 
well be nonnegligent at some distance from the threshold. 

8. The Bose statistics, the isospin symmetry, the C-invariance and the crossing, being 
the exact properties of the amplitude derived according to the rules of the quantum field 
theory, restrict the amplitude dependence on variables, making some parameters vanish and 
relating the rest parameters. 

9. As a result there might be about 20 parameters for the near-threshold description of 
all channels and, in particular, no more than 15 for the case of the Tr~p — > 7r~7r"'"n process. By 
the unitarity conditions which relate the nonvanishing 7r27r isospin amplitudes only with the 
P31 and Pii waves of elastic nN amplitudes, the above parameters must be approximately 
real. Only few of 20 imaginary parts of the polynomial background are expected to be of 
importance. 

10. The discussed structures must be combined with the explicit term describing the 
OPE amplitude and isobar contributions. This will provide about 20 additional parameters 
if D waves of the tttt scattering and all structures of the vertices ttttNA, tttt A inrNN^ are 
important. 
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Table 1: Resonances and particles from PDG-96 ||70[ responsible for pole contributions to 
the low energy ttN ttttN amplitude. 



11. The cuts due to thresholds of another inelastic processes like the Stt production, the 
rj production, etc. are conspired by nearby isobars. 

12. The analysis of the paper makes it evident that isobar resonances saturate the 
existing data on total cross sections below 1 GeV. Therefore, the imaginary part of the 7r27r 
amplitude might be described by the Breit-Wigner form of isobar contributions and the 
discussed above parameters of the imaginary background might appear to be negligible. 



2.3 Resonance Contributions to ttN ttttN Amplitude 

The 7r2'7r amplitude to be developed must cover the energy interval from the threshold to the 
isobar region. Therefore, the approach 2 mentioned in Introduction is the most suitable one 
for the discussed purpose since it admits the due account of isobar physics and can provide 
the maximal model independence of the obtained results. 

Let us first enumerate the resonance-exchange graphs which are proper to our reaction 
(it is suitable to use the two-particle channel for labelling). Three representatives SEttat, 
SE^TT, SEatat of inequivalent crossing-related single-pole graphs are shown in Fig. 1. Every 
four-particle vertex of the above graphs can be expanded in at most two different ways 
providing no more than three types of inequivalent double-pole graphs which we mark in 
DE7r7v,7rAf, ^'^nN,-Kn, DEata^.tttt (see Fig. 1). In what follows we shall use the self-explanatory 
notations like SE^Ar{A}, DENN,mr{^, p} when discussing contributions of classes of crossing- 
related graphs. It is evident that all distinct resonances and particles responsible for pole 
contributions to the amplitude enter already the single-pole scheme; the particles which 
might be relevant to the intermediate energy amplitude are listed in the Table |I|. There 
are enough representatives of the lowest spin-isospin states in the lists SE^^r, SEatat- In 
the case of the SE^ttv channel the list is far from being complete. In the present paper we 
limit ourselves by contributions from particles N, A, A^*, a, p, tt, uj and A; in fact, the 
interactions of a, u and A given below in the lists of lagrangians for the purpose of the 
forthcoming discussion were omitted from the actual analysis. 

Apart the OPE contribution our model is described by the effective interaction lagrangian 
which is used to construct the tree-level amplitude; all terms of the lagrangian are collected 
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Figure 1: Representatives of graphs of resonance contributions. 

in Tables |], |^. In general, we are trying to use the minimal derivative coupling especially 
when particles of the nontrivial spin (for example, in the pN A vertex) are being involved. In 
the contrary case the terms of too high order in the momentum are inevitable to be present 
in the amplitude. However, in the case of the three-particle vertices containing the pion, 
like irNN, tcNA, etc., the vertices are brought to the derivative-coupling form. Having in 
mind the complete equivalence of nonderivative and derivative couplings in such vertices (in 
the difference of the corresponding amplitudes the propagators always become contracted) 
we assume that only other types of vertices (for example, Air, irirNN, etc.) are responsible 
for the explicit breaking of Chiral Symmetry. 

Even in the brief summary of properties of the considered interactions we must point 
that the form of the four-particle vertices is chosen to represent all spin-isospin structures 
of the given vertex. This makes impossible to fix the lagrangian parameters from the decay 
data (for example, the known width for N^, ttuN provides only bounds for 4 parameters 
of the TTirNN^, lagrangian given in the Table 0). However, we prefer to avoid any doubt in 
respect to a possible model dependence of the results for vrvr-scattering lengths or, at least, 
to leave a chance to check the presence of such dependence. 

The central contribution to our amplitude is supposed to come from the OPE graph 
SEArAr(7r). It is parametrized basing on the cross-symmetric threshold expansion of the Air 
vertex elaborated in the paper . Apart the account of the imaginary part the expansion 



in the 0{k'^) order is equivalent to the form described in the papers |T^, |2T|. The results of 
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Graph 
[Vertex] 


Lagrangian 


[nnNN] 




[nnNN,] 


(7L^^.iV5,,iV,7r'^7r^ + gl^j,^^N6abN,d'^7r'^d,7r'' 

+gLNN,Nteabcr'7t.N4d'^7T''7T» - TT'^d'^Tr'] 
+(7^^^^^iV2e,,,r^[7^,7.]iV*9^7r«9'^7r^ + {H.C.} 


SE wiAl 

IrnrNA] 


NigLNAKc + ^aLNAFDi^^Aid^n^n' 
+N(ql maFL + 3q^ ^AFk^-f^-f^Afd^'n''d''n'' + |H.C.| : 

-^dbc = "^^dfec + ^dbTc - 5dcTb ; F}^^ = SdbTc + SdcU - 2/3 5bcTd 






[TTTTTTTr] 




SEjvArju;} 

iTnnruj] 

/I /I /I k-A^ 


a ,uj, d„7i"-dr.n''dfin''ie„hJe'"'"'^ 


SEAIAfMI 
'iViV L'' -^J 


jitttttA b ^^i" 

+gl^^^A^,d,7T'd,7T-d''n- + gt^^A'',d,7r%n-d'^7T- 


SE^^{o-} 




L J 




SE^^{p} 
[npNN] 


gl,NNNl,lul,N{d'^7^,p: - 5-7r,p^)/2 
+9lpNNNr''-f5Nd^Trbp'^ieabc 



Table 2: Lagrangians for 4-particle vertices of the SE graphs. 

the present work which will be discussed below confirm that the precision of the currently 
available experimental data does not allow to consider higher terms of the vrvr amplitude. 

We list in Tables |], ^ below the interactions most part of which had been used for 
the construction of the phenomenological amplitude entering the data fittings. Since all 
3-particle vertices of the SE series of graphs are present also in the DE graphs we do not 
describe them separately — see Table |^. In Table ^ the symbol V^tt stands for the vrvr 
interaction — this contribution is taken in the direct amplitude form. 
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Graph 
[Vertex] 


Lagrangian 


DE^iv,,riv{A^, N} 

[nNN] 




DE,;v,.;v{iV,iV,} 
[ttNN,] 


9.NN,NTa^^^l5N,^^'7^'' + {H.C.} 


DE^n,.n{N, A} 
[ttNA] 


(^^jvAlA^P^ftTVa/^Tr" + 7VP„\A»a^7r*) 








^.iv,A(A^P„;,7V,a'^7r'' + TV.pJ^A^a'^Tr'') 


[nN,A] 


DE7rAi-,7rJv{A, A} 

[ttAA] 


^,rAAA»7,75A„;,cA^9''7r'^ , 
^abc = |(— S^bcTa - S^acTft + S^afoTc — 5ieabc) 


DE^Af,^^{iV, a} 
[aNN] 


QaNNNNa 


[pNN] 


9pNnNJi,t^Np^^ + ^J^^iVa^,r"iVaV" 


[aNN,] 


^aArjv.A^A^*(J + {H.C.} 


DE^Ar_^jr{^*, P} 

[pNN,] 


^Mriv.A^7M^"^*P^ + ^,W^<^M.T«7V,aV'^ + {H.C.} 


DE^iV,7r7r{A,p} 

[pNA] 


^piVAA»P„biVp^^ + {H.C.} 


DEiV7V,7r7r{7r, CT } 

[ttttct] 




DEjV7V,7r7r{7r, p} 

[ttttp] 




[npLj] 




[iraA] 




DEjvAr,7r7r{^, P} 
[TipA] 

[ANN] 


5XiVAr^^7M75r'^iVAj^ + gl^^Na.^^.T'^Nd^^A'^ 



Table 3: Lagrangians for 3-particle vertices of the DE graphs. 
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2.4 Background Contribution to tiN ttttN Amplitude 

An important issue of the approach [|l^, ^ developed for the near-threshold energy region 
is given by the linear background terms presented in the form that respects the symmetries 
of strong interactions, namely, P, C, T, SUf{'^) and crossing. To modify this ingredient of 
the model we need first to discuss the role of these smooth terms in the amplitude. 

Initially, such terms were added to the OPE ones to stand for all other possible mecha- 
nisms of the TiN ttttN reaction. When taking into account only a part of contributions 
of particles listed in Table |l| one leaves enough room for the background terms standing for 
the rest resonances. 

Another reason for the presence of a background is connected with ambiguities which are 
specific to the off-shell interactions of high-spin particles (see, for example, the old discussion 
p6| , |5^ , |6^). These interactions result in polynomial terms in the considered amplitude (and 
polynomial terms in the 4-particle vertices). The overall contribution of this kind coming 
from all resonances is constrained by the asymptotic conditions for the entire amplitude. In 
principle, this makes necessary to use the consistent theory for particle propagators and all 
vertices (like irNN, ttNA, etc.). Leaving the parameters of the polynomial background free 
we are safe to use the propagators and the vertices in the form determined by the simplicity 
reasons and/or by the chiral-symmetry arguments. (In particular, we are using the simplest 
form 

S&{k) = l-SM^g'''^' + M^Y'Y + '2k^'k^ - M{k^'Y - Y' k^)] (27) 

for the nominator of the propagator of a (3^/2) particle.) 

Since we are regularizing the exchange graphs which have poles in the physical region 
by the iMT shift in the propagators the polynomial background must have both real and 
imaginary parts. 

The discussed nature of the background terms forces us to modify the model of the papers 
TP| , pT| in two directions. First, we add the second order terms in variables (|^ into the scalar 



structures Sa, Sb, Sc, Sd- 

= AuT^ + A^suu + A^Je + A^^iu^ + u^) + Aii,{e^ + e^) + A^jT{e + e) ■ (28) 
S^B^ = S^^\T,eu,eu,ee,ee); = ^^^(r, ez/, ez7, e^, e^) ; (29) 
5g) _ Aisiii^e - Ov) . (30) 

(This makes the dimensions of the tensor and the scalar structures of the decomposition (^) 
balanced.) 

Second, all 18 terms of the real background (11 of which are described in the paper 
pT| and the rest are given by eqs. (p8|)-(|30D) had been copied to provide the imaginary 
background of the amplitude. 

Let us now clarify the concept of contribution which we are widely using throughout the 
paper. From the point of view of Chiral Dynamics the usage of notions like the background 
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contribution (a part of which the so called contact terms are), the contribution of OPE, etc. 
is meaningless since the field redefinition does not allow a separate graph to be well-defined 
— see the relevant discussion in the book The absence of the common solution on the 
role of the higher-spin baryons in ChPT makes the above notions even more ambiguous. 
Nevertheless, such quantities as the residues of the poles and the on-shell parameters of V^t^ 
are well-defined — these very quantities are in the focus of our project. As for the off-mass- 
shell contributions, in particular, the OPE one, we take them as they are, the Lagrangian 
source and Feynman rules providing the model-dependent answer. The field redefinition 
then modifies, first, the polynomial background terms and, second, the parameters of SE 
graphs. It has been already pointed out that we are leaving the background parameters 
free (and keeping all spin-isospin structures of 4-particle vertices being represented) — this 
helps to avoid the dependence of the results on a particular model. 

To summarize, the considered amplitude contains 36 free parameters of the polynomial 
background, 4 free parameters of the real part of the OPE contribution and 5 formal param- 
eters of its imaginary part (their relations with parameters of the real part are described in 
the paper ||2^) and parameters from the lists of Table ^ and Table ^ discussed in the previous 
subsection. In the current paper we shall discuss the fittings performed with 21 parameters 
from the above lists coming from the exchanges SE^Ar(A^), SE^Ar(A), SE^Ar(A^*), SE^jr(p), 
DE^N,MN,N}, DE,jv,.^{Ar,A}, DE,;v,.;v{A^, iV*}, DE,;v,.^{Ar„ iVj, DE,;v,.^{iV„ A}, 
DE^7v,7rAf{A, A} only. For the purpose of the forthcoming discussion we need to enumer- 
ate the fitting parameters — this is done in Table ^ below. The background parameters 
Ai-Ais and their analogs iig-ise in the imaginary part of the amplitude have been already 
discussed in the previous subsection; these parameters had been being processed in fittings 
as they are — to save space we do not list them once more. Some factors were absorbed into 
the fitting parameters to improve running characteristics of the code. The list of expressions 
of the actual fitting parameters in terms of the constants of the interaction Lagrangian is 
presented in Table H. 
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Graph 


Fitting 
Parameter 


Expression 






J, AT J\T • Qo 




02 






■^o 












' 1 


yirpNN yTTTrp 




To 




DE \r W o\ 




V 

ilpNN I/TriViV i/TTTrp 




7-) 






7-) 


S'ttATA ■ 9mTNA 




n 


QnNA ■ 9-K%NA 




n. 


4 

Q-kNA • 9TrnNA 


i-'I^TTN,TTN\'-^- ^/ 




^ n r,r A Y'^ • n aa/Q 

^'jyynNAj i/TrAA/t' 


DE^AT^ArlA^, A| 

-^TTiV ,7riV L,-^ ' 7 ' — 'J 


^ 


1 Ott a 1 * Q'TT atm/o 


ctp f AT X 




S'TTAfAT* • 9\ttNN» 






9ttNN^ ■ 9-k-kNN^ 




p 


^ ^3 

'^9t:NN; ■ 9TnyNN, 




/?i 


C/„ \ ' 1- • c; ^ ,, , . 

y TTiV i\ * y ttttA a » 


DEjrAr,7rAf{^*, 


i?5 


9i^N»N» ■ {9'nNN,Y/'& 


DE^AT^^ArjA^*, A} 




9nNA ■ 9itN,A ■ 9TTNN,/'i 


DE,iv,.iv{A^,Ar,} 




9itNN ■ {9-nNN,Y/^ 


SE,^{iV} 




9ttNN ■ 91ttNN 






2 

9itNN ■ 9-kttNN 






'^9-kNN ■ 9I-KNN 






4 

9ttNN • 9TTnNN 


DE,^,,^{iV, iV} 




{9nNNy/& 



Table 4: Expression of fitting parameters of the phenomenological amphtude in terms of 
the Lagrangian parameters. 
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3 Experimental Base 



The experimental base for fitting the amphtude parameters was built of the total cross sec- 
tions in all five channels of the considered reaction in the energy region PLab < 500 MeV/c 
and the distributions measured in the hydrogen bubble-chamber experiments at these en- 
ergies pS, H6|, |T^, ITT], l3S]. The bubble-chamber data were preferred since they satisfy the 



condition of the coverage of the reaction phase space in the most complete way (the only 
restrictions are Pproton > 120 MeV/c and > 30 MeV/c which correspond to 5 mm of the 
flight distance). The losses occupy only a small part of the phase space (namely, < 2 % in 
the case of {n^n^p} and < 3 % in the case of {Tr^Tr+ri} at Pbeam = 400 MeV/c) and might 
be easily taken into account during the calculations of theoretical distributions. Besides, the 
systematic errors of the bubble-chamber experiments are also minimal. 

In the course of our fittings we faced the following problem. There are three works in 
the {iT^TT^n} channel at close beam momenta: 

1. Kirz (PLab = 485 MeV/c) ||, 

2. Blokhintseva (PLab = 457 MeV/c) [||], 

3. Saxon (PLab = 456 MeV/c) |7ll. 

It was found that results of these works are incompatible. If one tries to describe only 
distributions from these works without appealing to any other data (total cross sections or 
other channels' data) then the averaged P^r bin (= x'^p.b.) distributions by Kirz and 
Blokhintseva separately and for their sum as well is less than 1. Adding the distributions by 
Saxon we obtained x'^pb > 1-5 for every combinations and even x'^pb — 2-1 for the Saxon's 
distributions themselves. 

We should note here, that we have the results of the experiment by Blokhintseva as the 
collection of events. Hence, we can build any distribution and we did so for all kinds of 
distributions published in the Saxon's paper. Along with other spectra such distributions 
also represent the Blokhintseva work ]18|] in our fittings. 

The experiment by Saxon contains more events than the ones by Kirz and Blokhintseva 
together, so, at first glance, one should prefer to choose his results. However, the processing 
of films had been performed by the Saxon's group in a quite specific (nontraditional) way. In 
such films the elastic {vr~p} reaction events look very much like the events of the considered 
reaction, their total number being ten times greater. To economize manual measurements the 
group had estimated visually the density of the positive tracks and, basing on the estimate, 
had selected the events with vr"*". 

However, besides the velocity of the charged particle the density of a track in the bubble 
chamber depends upon much more another factors like the moment of particle flight relative 
to the moment of liquid expansion start, the moment of snapshot relative to the moment 
of particle flight, the liquid superheating degree, etc. This dependence is extremely strong. 
The above parameters undergo stabilization but the latter never becomes ideal. The rest 
fluctuations can not prevent the ionization measurements of tracks but for the reliability of 
the determination of the particle velocity in the every snapshot the control measurements of 
the bubble density should be performed for the sample track. In the visual estimating the 
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probability of an error is big especially in the case of the so large number of the processed 
tracks. 

One can estimate the systematic error of the Saxon's work in the following way. The 
elastic events accidentally selected into the list of the considered reaction anyway will be 
rejected after measurements by the reaction fit results. The inelastic ones for which the 
positive track was erroneously identified as the proton one are lost forever. In such a way 
the statistics become more poor, the parts of the phase space for which the vr"'" velocities 
are small suffering the most significant losses. The measured value of the total cross section 
is given in the Saxon's work as 1 mbn; since the results of all isotopic analyses of all set of 
total cross sections data provide the value of 1.4 mbn one can deduce the 30% level of losses 
of events in the "economizing" routine of the reaction analysis. 

So large value of systematic error forced us to withdraw the distributions of the paper 
| 7T| from our fittings. 

Thus, the distributions liable to treatment belong to three channels of the considered 
reaction: {7r~7r"'"?7,}, {7i~7i^p}, {vr+TT+n}. The data on distributions of works ^ 



were taken from the journal publications. The major part of distributions of the work 
ITSl has never been published. These are the distributions in the following variables: the 
squares of invariant masses of all pairs of final particles W^— Wj^-n, Wj^+n, the invariant 
variables r, uji, uj, On, 6j introduced in [|19], the cosine of the angle between the planes 
determined by CMS momenta of a) the beam and the recoil neutron; b) vr"*" and vr" (cos6'^7r); 
the cosines of CMS angles of final particles with the beam cos^t^-, cos^^+, cos6'„, the angle 
D(f) = 0~ — 0+ of the planes determined by CMS momenta of a) the beam and the tt" meson; 
b) the beam and the vr"*" meson, the azimuth angle (phe&m and the cosine of the polar angle of 
the beam cos^beam in the Saxon reference frame in which the x axis is the direction of the 
neutron momentum, the z axis being along the vector product of the neutron momentum 
with the momentum of 7r+ meson, [Pn,P7r+]- 

The Jones paper (the {vr^vr+n} and {n^TT^p} channels at P^ab = 415 MeV/c) pro- 
vides distributions in the squares of invariant masses of all pairs of final particles (namely, 
{vr-vr+n}: W^-^+, W^-n, W^+n, {n-n'^p}: W^-^o, W^-p, W^Op). 

The Kirz paper (the {7r~7r+n} channel at pLab = 485 MeV/c) gives distributions in 
the invariant masses 1^^-^+, W-^-n, W^7r+n of all pairs of final particles and in the cosine of 
angles of all final particles with the beam in CMS co^O^^-^ cos^^+, cos6'„. 

Another paper by Kirz (the {7r+7r+n} channel at pLab = 477 MeV/c) contains dis- 
tributions in the CMS angles of final particles 6^^+, 6'„, the tt"*" CMS energy T^+cms? the 

momentum of vr"'" in the (Tr"*", vr"*") system and in the nonrelativistic momentum transfer 
p 

-' transf- 

One should note that some of these distributions strongly differ from the distribution 
provided by the empty phase space while there are some with the insignificant difference. 
This becomes especially clear in terms of (normalized) quasi-amplitudes — see pictures in 
Figs. 3-8 where almost all distributions are reproduced. 

We are using the total cross sections data §, §, §, [|], §, 0], [0, [0, 0, 0, [^, 

H, m, n, Q, n, n, m, h, q, ii, h, 0, o, g, Ei> El- 94 points 



were selected on the grounds of the compatibility analysis of the paper from the total 
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list of 105 experimental points. 

In some variants of fittings (namely, when the theoretical amplitude was not completely 
real — see the next section where all variants are described in more details) there were used 
two experimental data points more. These points were fixing the phases of two (nonvanishing 
at the threshold) isotopic amplitudes by the known values of the elastic P31 (~ —4°) and 
Pii (~ 2°) phases in accordance with the final-state interaction theorem. 



4 Data Analysis 

The analysis of the data described in the previous section required some preliminary steps 
which are described below. After then we discuss the main results in the subsect. 4.2. 



4.1 Principal Steps of Analysis 

To perform the analysis the following steps had been done: 

1. The contributions to the amplitude of every parameter of our model (discussed in 
the sect. 2) had been calculated analytically according to the Feynman rules of the tree 
approximation. (The iMT shifts regularizing the terms with poles located in the physical 
region were made also for all cross terms to ensure correct crossing properties of the entire 
amplitude.) The contributions generally come to all 16 scalar-isoscalar form factors defined 
in eqs. (H), (HI); the corresponding expressions had been obtained and transformed to the 



FORTRAN code with the help of the REDUCE package for analytic calculations in high 
energy physics. 

2. The second step was aimed to simplify the data fitting routine and save a lot of time 
during thousands of fitting iterations. For a given experimental point (a) the theoretical 
quantity confronting the experimental value a^^^ is the integral over the reaction phase 
space (or over its part in the case of the distribution data) of the squared modulus of the 
amplitude (p^Ol). 

Since the set of all parameters {A,y} including the formal ones enters our amplitude 
linearly and the conditions for the formal parameters do not contain kinematics we can 
present the theoretical quantity ([T2|) in the form 

fl,U 

Then we can build for every experimental point (a) the correlator C^^-^ performing all phase 
space integrations ones and forever. These calculations had been done for every bin of all 
distributions with the use of the standard high energy physics package based on the Monte 
Carlo integration. Since all distributions in question belong to the beam momenta 415, 460, 
477, 485 MeV/c few Monte Carlo runs were necessary (and sufficient) to provide the calcu- 
lations. In the case of the total cross sections the correlator matrix C(a) for an experimental 
point (a) was being recovered during the fittings' run-time by the fast interpolation from 13 
fulcrum matrices in every channel precalculated at beam momenta 280, 285, 290, 300, 325, 
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350, 375, 400, 450, 500, 550, 600 and 650 MeV/c. The calculations of the fulcrum matrices 
had been performed with the fast and efficient program of Gauss integration described in 
the paper |2§. 

3. The content of the third step was the data fittings themselves. Here, we point only 
the general features and specifics of the approach; more details will be revealed along with 
the discussion of results. 

The fittings had been performed by minimizing the value of defined as 

(a) ) 

in the case of the total cross section data. However, to avoid the artificial increase of the sta- 
tistical weight of the total-cross-section data point for which the numerous distributions also 
were undergoing fittings we were using the distributions' data cr^afn which were normalized 
to 1 instead of the total cross section. 

The use of the precalculated correlator matrices C*(q,) and the simplicity of the expression 
(0), being calculated in the course of the iterative minimization of x^? made it possible to 
perform thousands of fitting runs each of which executing hundreds or even few thousands 
of iterations. The same specifics of the approach provided the excellent flexibility in respect 
to variation of both the set of fitting parameters and the set of experimental points. The 
numerous variants of fittings will be discussed in the next subsection. 



4.2 Major Results of Analysis 

Already test runs of fitting the distribution data had shown that the simple model of the 
paper is unlikely to be capable to provide a satisfactory description. Therefore, the first 
question we were trying to find the answer to was if there exists a relatively simple model 
consistent with the selected data base (411 experimental points). The answer was "no" and 
let us now discuss why. 

We had been grouping together parameters related to one or another mechanism of the 
considered reaction since it ought to be a hard task to test ~ 2^^ variants of all combinations 
of parameters. There were 7 such groups in total, the largest ones consisting of the sets of 
the background parameters. Below we are using the following symbolic notations for these 
groups: 

"b" — parameters of the real background; 
"i" — parameters of the imaginary background; 
"o" — parameters of the OPE contribution; 
"r" — 3 parameters related to the p meson; 

"N" — parameters of the nucleon contributions (namely, SE^7v{A^} and DEt,j^^t,j^{N, N} 
ones) ; 

"D" — parameters of the A contributions (namely, SEt^tvIA}, DE^atttAtIA^, A} and 
DE^Ar,^Ar{A, A} ones); 

"R" — parameters of the A^* contributions (namely, SE^ArjA'*}, DE.„n^t,n{N, N^,}, 
DE^ N,nN{N^,N^} and DE^Ar,^Ar{iV=^, A} ones); 
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For all possible combinations of the listed groups there had been performed the data 
fittings with at least 50 random starts. In general, no new solutions had been being revealed 
after 20-30 starts. In the cases when the number of distinct solutions was greater than usual 
we had been continuing the search with the increased number of random starts to 100 and 
more. With one exception this had not provided new solutions (apart the solutions with 
incomparably large x^)- Only in the "DRNbior" variant when all groups had been being 
involved acceptable distinct solutions had been found after 100 starts. The total list of this 
variant contains 12 solutions, not all of them are shown in the discussed Table. (It should be 
noted that the latter variant is the most difficult in fitting: the convergence is slow because 
of huge correlations between parameters, the parameter errors in the gained solutions being 
too large.) Therefore, we consider the probability of existence of a missed solution to be 
negligible. 

The values of of the best solutions are listed in the Table ^. A part of variants is 
withdrawn there to make the table more compact. The baryon-exchange mechanisms are 
combined with the sets of parameters "b", "i", "o" and "r". For example, the box in the 
"DR" row and "bio" column corresponds to "DRbio" variant, etc.; in the column (row) with 
the empty title we collect the results corresponding to the "pure" mechanism. Initially, the 
"r" set was considered as a perturbation to the basic amplitude since the relatively narrow 
p resonance might get into the reaction phase space at PLab ~ 800 MeV/c only. (The use of 
the higher polynomial background in such a role requires to make calculations with about 
a hundred of the next order terms with free parameters or to find the reasons for rejecting 
the most of them.) However, the presented values of can not provide the inference that 
the "r" mechanism is unimportant (unless all baryon-exchange mechanisms are committed 
to action). 

It is the examination of the Table ^ which leads to the following conclusions: 

1 . There is no particle for which a simple exchange mechanism is capable to describe the 
data. 

2. Among the double-particle exchanges the participation of OPE does not look advan- 
tageous. (Even the dummy "r" mechanism looks more preferable sometimes.) 

3. The leading order ChPT fails to describe the data. (Indeed, all contact terms of 
Chiral Dynamics and tree-level graphs are contained in the variant "Nbo" ; the result of the 
latter even for parameters being free looks depressing.) 

4. The most significant improvement of is achieved when the imaginary background 
"i" is being added; the participation of A and iV* exchanges with the strong imaginary part 
of contributions is leading almost to the same effect. 

To make the above general conclusions the knowledge of the absolute values was suffi- 
cient, especially, since the number of experimental points considerably exceeds the (varying) 
number of free parameters. For a more subtle deduction one needs an information on Xbf 
(i.e. per degree of freedom). However, before calculating Xbf first ought to look if 
there are undetermined and inessential parameters in the solution. 
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Table 5: List of in distinct solutions obtained from random start (50 tests) with 
A^Exp = 411 experimental points. (In the truncated variant "DRNbior" there are 12 solutions 
in total obtained in fa 500 runs.) 



Now we must note the following feature of the discussed results. The numerous solutions 
in some variants reflect the existence of large correlations of parameters in the variant in 
question. The estimating of parameter errors by the fitting routine is not then precise enough 
and our simple algorithm for comparing solutions marks as different the solutions which are, 
in principle, identical. 

To make the situation clear and, what is more important, to reduce the errors of parame- 
ters we had performed the series of additional runs eliminating one by one those parameters 
which relative errors had been being greater than 1 and fitting the data by the rest pa- 
rameters (while the value of Xdf been improving). This routine had been applied to 
all solutions. Unfortunately, there is no enough room to display the results. The Table ^ 
presents the information on Xbf some selected variants, the number of parameters in 
effect A^p and the independent lowest vrvr scattering lengths. The results of eliminating pa- 
rameters had no influence on the above conclusions 1.-4. whereas the number of parameters 
in effect had been considerably reduced. The new inferences derived from the Table ^ read: 

5. The considered data base requires a complicated model for its description. Below the 
level of Xdf < ^-^^ absolute minimum of the parameter number is 22 (variant "Nbi"); 
for the most part of the acceptable solutions the number is greater than 30. 

6. Assuming the consent that the sequence of signs +, — , + of the tttt scattering lengths 

physical, one finds that a half of all acceptable solutions with the OPE 
contribution falls to the unphysical sector. (The interpretation of the latter phenomenon 
in terms of threshold values of amplitudes will be discussed in the next section. Here, we 
simply note that the physical sequence is the one consistent with the predictions of Chiral 
Dynamics.) 

7. The results of fittings can not improve the precision of determinations of tttt scattering 
lengths. For example, the region 0.06 < Oq^" < 0.19 might be derived from the Table ^ as 
the preferable one; however, even the current experimental value ||3^ ag^^ = 0.26 can not 
be rejected on the ground of the Xdf criterion. 



8. The predictions of ChPT in the next-to-leading order |^ Oq^" = 0.20, Oq^^ = —0.041, 
a{^^ = 0.036, 02^° = 0.002 are compatible with the data base; the separate fittings with 
the OPE parameters being kept fixed by the above values of scattering lengths, resulted in 
solutions with the values of Xdf ^^ry close to the best ones; for example, in the "DRNbior" 
variant we get Xdf range 1.175 — 1.200. 

9. The formal parameters of the imaginary part of the OPE contribution are insignificant 
at the explored energy region (this is the result of the separate investigations; eliminating 
these parameters we get some shifts in the parameters of the imaginary background and in 
the parameters of A^* and A isobars, the real parameters of OPE remaining the same). 

The discussed above routine of eliminating parameters provided further support of the 
point 5.: no one of the considered mechanisms had been rejected as a whole in the course 
of improving the Xdf value. The listings of this routine contain also a large volume of 
information about the significance of a given parameter in terms of the "statistical" frequency 
of its participation. For example, in all variants the parameters A2, A3 and A4 of the linear 
background of ref. and the parameters and Ai^ of eqs. (p8D, (p9|) were found to be 
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bi 




bio 










bir 




bior 












9 




9 

-Vdf 




1=0 

"o 






2 

Adf 


Np 


2 

-Vdf 


Np 


"o 








1.94 


26 


1.56 


30 


.04 


-.151 


-.016 


1.58 


31 


1.46 


32 


.009 


.15 


-.040 




1.95 


26 


1.57 


29 


.03 


-.155 


-.011 


1.58 
1.64 
1.66 


32 
31 
31 


1.46 
1.50 
1.51 
1.53 


32 
30 
33 
33 


.014 
.136 
.044 
.045 


-.16 
-.12 
-.16 
-.20 


.020 
.048 
.048 
.039 


D 


1.45 


33 


1.37 


24 


.16 


-.084 


.040 


1.40 


30 


1.31 


30 


.14 


-.14 


.056 




1.45 


25 


1.40 


30 


.11 


-.187 


.064 


1.42 


38 














1.45 


27 












1.42 
1.44 


36 
31 












R 


1.35 


32 


1.33 


32 


.11 


-.004 


-.026 


1.25 


37 


1.22 


36 


.038 


.015 


.000 




1.36 


31 


1.33 


33 


.12 


.002 


-.025 


1.38 


37 


1.35 


39 


.073 


.0002 


.008 








1.51 


35 


.18 


-.138 


-.012 






1.40 


34 


.170 


-.034 


.028 


N 


1.42 


23 


1.39 


27 


.07 


.004 


-.014 


1.36 


33 


1.32 


30 


.23 


.052 


.011 




1.42 


22 


1.39 


28 


.07 


.004 


-.014 


1.37 


27 


1.32 


30 


.23 


-.012 


.018 








1.40 


33 


.16 


.026 


-.020 


1.37 


31 


1.32 


34 


.14 


-.012 


.019 








1.41 


28 


.16 


.023 


-.020 


1.37 
1.37 


28 
27 


1.33 

1.33 


29 

26 


.23 

.10 


.029 

-.014 


.017 

.015 


D 


1.25 


34 


1.25 


30 


.00 


.000 


.000 


1.20 


36 


1.21 


41 


.13 


-.026 


.046 


R 


1.27 


33 


1.26 


35 


.19 


-.058 


.038 


1.26 


40 


1.24 


46 


.12 


-.105 


.057 




1.29 


35 


1.26 


35 


.10 


-.048 


.033 


1.26 


39 














1.33 


38 


1.31 

2.01 


43 
40 


.10 

.80 


-.016 
.21 


.026 
.100 


1.28 


38 












D 


1.31 


31 


1.29 


32 


.04 


.025 


.005 


1.27 


34 


1.19 


35 


.11 


.000 


.045 


N 


1.31 
1.36 


32 
26 


1.30 


32 


.15 


-.022 


.023 


1.29 
1.33 


34 
33 


1.24 


31 


.17 


-.104 


.048 


R 


1.24 


29 


1.23 


36 


.22 


.037 


-.035 


1.19 


35 


1.15 


43 


.14 


.042 


-.008 


N 


1.34 

1.39 


39 

36 


1.31 

1.32 


36 

36 


.09 

.20 


.009 
-.047 


-.015 
.0002 


1.32 


45 


1.29 


40 


.18 


-.048 


.005 


D 


1.18 


41 


1.17 


39 


.05 


-.027 


.090 


1.16 


44 


1.14 


47 


.15 


.041 


-.026 


R 


1.20 


36 


1.19 


35 


.02 


.017 


.005 


1.18 


41 


1.15 


46 


.07 


.025 


.028 


N 


1.21 


36 


1.22 


38 


.27 


.005 


.028 


1.19 


40 


1.16 


42 


.07 


-.056 


.045 




1.22 


37 


1.22 


38 


.05 


-.029 


.015 


1.21 


36 


1.20 

1.20 
1.21 


45 
41 
39 


.07 

.17 
.19 


-.076 

-.053 
-.059 


.047 

.054 
.053 



Table 6: List of Xbf distinct solutions obtained after deleting undetermined parameters; 
Nexp = 411 experimental points. (To save the space only the best unphysical solutions of 
the "DRNbior" variant are shown.) 
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necessary whereas the parameter An entering the tensor structure of the amphtude D had 
been always ignored. In respect to the OPE parameters the content of the Tables | already 
shows that the parameters are not of the top significance since there are many acceptable 
solutions without the OPE contribution at all. The parameter gi of the paper |^T| is found 
to be relatively important, the Z)-waves parameters g2, gs being much less necessary. 

For illustrations we have chosen the solution from the "DRNbior" variant with Xbf = 1-16 
(see Table The data on total cross sections and the theoretical curves in terms of the 
quasi-amplitude ([15|) are drawn in Fig. 2. The most intriguing feature of the discussed 
curves is expressed by the practical coincidence of the theoretical results for the {—Op} and 
{+ Op} channels — this is clearly seen from the separate picture in Fig 2 which we draw for 
the combined data. The picture makes it obvious that the discussed phenomenon is strongly 
implied by the experimental data. This means that the isotopic amplitude D, being the only 
origin of the difference of the theoretical cross sections of these channels (see eqs. @), must 
vanish indeed. 

At Figs. 3-8 one can find the theoretical curves of the same solution for the distribution 
data discussed in the sect. 3 (in terms of the normalized quasi-amplitude (|T3p). A part 
of the data did not enter the fittings — in such cases the curves actually represent our 
predictions. The normalized quasi-amplitude measures the deviation from the empty-phase- 
space pattern. In its terms the behavior of the experimental data themselves appears to be 
quite different. 

For example, almost all angular spectra look flat. The explanation is simple: unlike 
the case of elastic reactions 2^2 the ID distributions are given by the remaining (3- 
dimensional) phase space integrals in our case of the 2^3 process. The averaging over 
polarizations together with this averaging over the reaction phase space make the angular 
dependence so weak. However, this can not be true for sections of the phase space and/or 
for higher dimensional distributions. 

In contrast, the distributions in the invariant variables ([^) are nontrivial; this proves the 
choice of variables of the paper to be the characteristic one for the dynamics of the 
considered reaction. 

The results of data analyses described in the current section make it possible to use 
the obtained amplitude for modehng two other approaches discussed in the Introduction, 
namely, the one by Olsson and Turner and the Chew-Low extrapolation. What we learn 
from these tests is discussed in the sect. 5. and 6. 
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5 Modeling the Olsson— Turner approach 



The original idea of the Olsson-Turner approach was to relate the threshold values of the 
7r27r amplitude with the vrvr scattering lengths. Therefore, there are two principal steps 
in the discussed approach: 1) determination of the threshold limits of amplitudes of the 
independent isospin channels; 2) calculation of vrvr scattering lengths with the account of 
other contributions to the above threshold limits. 

5.1 Threshold AmpUtudes 

Let us discuss the first step. 



5.1.1 Database 

We already know that the threshold amplitudes of 5 experimentally observable channels 
are not independent. To provide the experimental test of the relations (|2^) the threshold 
limits \M^\ must be determined from fittings of the data on total cross sections. The 
results of the similar fittings were already reported starting from publication (see also 
p7|). The most important conclusion of the work |]2^ is that the data in the region of 



-^Lab < 400 MeV/c admit amplitudes which are linear in the CMS energy. Here, the data up 
to Phab < 500 MeV/c will be exposed to linear fittings. Moreover, any preliminary selection 
will be excluded. Indeed, it is the unmotivated preference of one or another set of data which 
is the reason for contradictory results. 

The entire database contains along with the old data (and very old ones) the relatively 
new results. It should be noted, that before the OMICRON measurements [^-[^ only the 
{ — \-n] channel allowed to obtain the definite results of the linear fit; the above mentioned 
OMICRON data provided the possibility to carry out the procedure also for the {—Op} and 
{+ + n} channels. The authors of the paper [2^| took advantage of the precise data on 



the cTjoon} (very close to threshold) and provided the simultaneous hnear fit of all channels. 

The recent experimental information |^ for the first time makes possible to determine 
the threshold limit of cross section of the channel {+0p} and to test the prediction of the 
relation 

M°+op} = Mf.o,! . (33) 

(In fact, the approximate equality of cross sections of these two channels is observed along 
all energy interval considered.) 

The separate attention we pay to the {+ + n} channel. There are two sets of the near 
threshold data in the channel namely, |^ and which are in certain disagreement. 

They lead to different threshold limits |M|^^^||. 



5.1.2 Procedure 

Therefore, to provide the definite conclusion four basic variants of fitting were used. Their 
symbolic notations, entering the summary Table |^, are: ALL — variant with all the data 
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being included; OMI — in this variant the data are omitted; TRI — variant with the 



exclusion of the data [Q; X — the variant in which both sets |^ and are excluded 



The notation IND is used for the case when the data in each channel are fitted separately 
(the channel {+ + n} was treated in accordance with the selection adopted in ALL, OMI, 
TRI and X variants). This helps to understand the trend dictated by data of individual 
channel in the simultaneous fit. 

There were two types of fits in all four basic variants. First, the independent threshold 
values and are treated as being real (or, the same, as obeying the trivial relative phase). 
We use 7 parameters for simultaneous linear fittings of 5 channels: 2 independent parameters 
for threshold values (namely, and 5°) and 5 independent slopes in the invariant kinetic 
energy 

Tk = - To . (34) 



(In the course of calculations the latter quantity was taken to be Tk = y{p + kiY — {mj + 
1^2 + /^s) since isospin splitting in the particle masses can not be processed as a correction 
because of the nonanalytic dependence of the near-threshold phase space on masses — the 
demonstration can be found in [0.) 

Second, we add one parameter to describe the relative phase of the complex quantities 
A^ and 5°. 

There is the strong motivation to consider the parameters A^ and real. The estimate 
of the imaginary part of the 7r27r amplitude by the dispersion analysis of the paper ^ allows 
to consider the amplitude to be approximately real up to the energies pLab = 0.50GeV/c. 
From the general point of view of unitarity relations the value of the imaginary part collects 
contributions which are due to the following: 

a. Three-particle intermediate states. This includes contributions of quasi-two-particle 
states when the third particle is present as a spectator (namely, for the vrvr state and tt 
for the ttN state; the former configuration is characterized by the large imaginary part of 
isospin-zero vrvr amplitude). The three-particle phase space makes this contribution vanish 
at threshold. (In the paper with the use of (1/m) expansion it was already verified that 
TTTT loops did not contribute to the threshold 7r27r amplitude.) 

b. Two-particle intermediate states. In this case only the vrA^ system is allowed. The 
angular momentum conservation forces these particles to be in the P wave when the final 
particles at threshold are in the S wave according to Quantum Mechanics — otherwise the 
P-parity properties of the vrA^ system mismatch those of the ttttN final state. The known 
phases of the Pu and P31 waves of the vrA^-elastic amplitudes approve the neglect of the 
imaginary part of the tt2tt amplitude at the threshold. 

c. Single-particle intermediate states. The only possible one at the threshold of the 7r27r 
reaction is the A isobar. It appears in the P wave of initial particles and its decay into the 
threshold configuration of the final vrTriV state has negligible width ||7CI|| . 

Therefore, there should be no expectations for finding a physically meaningful imaginary 
part (or a nontrivial relative phase of A^ and P°) at the 7r27r threshold. Hence, the threshold 
identity ( pS] ) discussed in the subsect. 2.1. must be practically exact. (The identities ( pT]) are 



the exact consequences of the isotopic invariance irrespective of the value of the imaginary 
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part.) 



5.1.3 Results 

The experimental data are shown in Fig. 9 where the values of quasi-amplitudes {Mx) 
are plotted versus the invariant kinetic energy Tk and the threshold amplitudes \Mx\ are 
collected in the Table |^. 

This Table is organized in the following way. The first column contains the symbolic name 
of the selection used for fitting. Five subsequent columns refer to the channel described in 
the column header. The last column provides the characteristics of the fit — Xdf- To make 
the comparison easy the additional row IND contains the results of individual fits of every 
channel (in the case of the {+ + 72} channel three additional fits with the exclusion of the 



controversial data sets, namely ||72[ and [|^, were performed). The last two rows present 
the predictions of eqs. (|25|) and (^) based on the data of individual fits of the channels 
{- + {OOn} and {-Op}. 

The bold-face numbers in the boxes ALL, TRI, OMI, X display the results of the solutions 
in the main fit with 7 = 2 + 5 parameters. Every time there was also present another solution 
of the fit (the slanted numbers in the Table |^ are used for resulting values). It originates from 
the sign ambiguity which was resolved in part by means of the phenomenological analysis 
of the subsect. 2.1.3. These auxiliary solutions choose the relation ( p^ to be true one 
and it was called unphysical since it provides the equal signs of isospin-0 and isospin-2 vrvr 
scattering lengths. 

The third line in every box of variants ALL, TRI, OMI, X represents results of the fit 
with the additional parameter describing the relative phase of amplitudes and when 

and are considered as complex numbers. The resulting values of the imaginary parts 
of the amplitudes M|_^^| and M|qq^| are given in brackets in the fourth line (the overall 
phase ambiguity is resolved be assuming the amplitude B^ to be real). 



5.1.4 Discussion 



First of all, let us compare the results of the linear fit discussed here with the results of the 
previous section — it is sufficient to examine the pictures on Figs. 2, 9. One can deduce 
that the threshold limits provided by the linear fit are good for all channels but the { — Vn] 
one. Because of the perceptible curvature revealed by solutions in the latter case the linear 
fit generally underestimates the value in question. 

The Table |^ does not allow to make an unambiguous conclusion on the ground of the 
quality of the fit only. According to the criterion all solutions are practically on equal 
footing. Even the unphysical ones can not be formally rejected. 

The relatively low values of of the unphysical solution are, first of all, due to rather 
large overall uncertainties of experimental data. The more discouraging reason is the exis- 
tence of arrays of data (almost in every channel) supporting the discussed solution. Only 
the data of the {0 n} channel as a whole reject this solution. The examination of the defi- 
nitions ( p^D displays that at the value of B^ being fixed by the data of {— Op} and {+ + n} 
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Fit 


{- + n} 
\A^ + B°\/V2 

\ ' \ / V " 


{ n} 

l^°l/2 
1 1/ 


{-Op} 

1 1 / 


{+0p} 

1 1 / 


{+ + n} 


Xdf 


ALL 


398 ± 18 

524 ± 24 
401 ± 150 


413 ± 12 

252 ± 12 
409 ± 101 

/ c:9 _L 07 
(JZ zc / J 


132 ± 6 

J18± 5 
132 ± 6 


132 ± 6 

118 ± 5 
132 ± 6 


264 ± 11 

236 ± 11 
264 ± 11 


0.99 

1.56 
1.00 


TRI 


404 ± 18 

518 ± 24 
406 ± 26 


406 ± 12 

260 ± 12 
406 ± 18 


121 ± 6 

106 ± 6 
119 ± 6 


121 ± 6 

106 ± 6 
119 ± 6 


241 ± 11 

212 ± 11 
238 ± 12 


0.78 

1.33 
0.77 


OMI 


372 ± 23 

533 ± 32 
401 ± 60 
( 181 ± 23 ) 


447 ± 13 

241 ± 13 
410 ± 54 
( 128 ± 16 ) 


184 ± 10 

136 ± 10 
197 ± 11 


184 ± 10 

136 ± 10 
197 ± 11 


368 ± 20 

272 ± 20 

395 ± 22 


0.83 

1.60 
0.80 


X 


387 ± 26 

506 ± 36 
402 ± 98 
( 154 ± 40 ) 


427 ± 14 

275 ± 14 
410 ± 76 
( 109 ± 28 ) 


154 ± 12 

83 ± 12 
164 ± 14 


154 ± 12 

83 ± 12 
164 ± 14 


307 ± 24 

166 ± 24 
329 ± 29 


0.76 

1.36 
0.76 


IND 


401 ± 19 


409 ± 18 


155 ± 25 


97 ± 51 


263 ± 11 
237 ± 12 
429 ± 26 
358 ± 37 


ALL 
TRI 
OMI 
X 








125 ± 22 


125 ± 22 


250 ± 44 




Eq.dH) 








155 ± 25 


310 ± 50 





Table 7: Threshold values provided by data fittings in the variants: ALL — all data; TRI — 
without OMI — without X — excluding |72]. Bold-face numbers correspond 
to physical solutions of real fits with 7 = 5 + 2 parameters, slanted ones — to unphysical 
solutions. The roman font is used for fits with complex amplitudes. The numbers in brackets 
in the columns of the { — hn} and {OOn} channels display the imaginary part of the resulting 
amplitude (provided is real). 
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channels at the scale no less than 200[GeV]~^ the threshold limits of { — h n} and {OOn} 
channels can not be in balance provided and amplitudes are of equal signs. While 
the closest to threshold points of the OMICRON experiment well agree with the high 
value of |M|__(^^|| the data of the {OOn} channel |^| leave no room for describing the cross 
sections of this channel at so small level (see Table |^). 

We end up our discussion of the unphysical solution by noting that there is a striking 
feature of the entire database: the distribution data, being added to the total cross section 
ones, did not help to make the distinction between these two kinds of solutions. Moreover, 
the absolutely best solutions had been found to be the unphysical ones — see Table 

In all variants the presence of the imaginary part improves the fit. In the variants 
ALL and TRI the imaginary part Im remains consistent with zero: 104 ±174 — ALL 
and 10^^ ± 10^ — TRI (the amplitude is considered to be real). In the rest variants 
the imaginary part is found to be unreasonably large: 255 ±32 — OMI , 218 ±57 — X. 
Nevertheless, it can not help to get rid of the contradiction of data in the {—Op} and 
{± ± n} channels, especially in the OMI variant of the individual fit. 

The examination of the threshold identities (0), ( P5| ) clearly displays that the contradic- 
tion is between the OMICRON data in {— p} and {±±n} channels themselves. Therefore, 
the analysis of the Table |^ provides the main conclusions that the most consistent from the 
point of view of threshold identities, the stability against the addition of the complex phase 
and a reasonable distinction from the unphysical solutions are the threshold amplitudes given 
in the TRI and ALL variants of the above Table. 



5.2 Account of NonOPE Contributions 



The discussed above threshold amplitudes are the keystone of the Olsson-Turner approach. 
(Therefore, it was the inconsistent input of the OMICRON analysis |^l|-[0 that provided a 



controversy in the derived values of vrvr scattering lengths.) The original formulae |^ relating 
the 7r27r threshold amplitudes with the vrvr scattering lengths were based — in modern terms 
— on the leading order Lagrangian of ChPT. 

Now the improved formulae take into account the next-to-leading order terms of the 
Chiral vrA^ Lagrangian |lT|. The series of papers [jl2|, ^ |l3l deals with various schemes 
of accounting the above terms in the framework of ChPT and/or Heavy Baryon ChPT 
and gives the predictions for the nonrelativistic quantities Di and D2 which up to a factor 
coincide with the discussed above threshold amplitudes and B^: 



Do 



2m 
m+Eo 



B' 



'4mv^l A° 



(35) 



The difference in the predicted values (which are quoted in the Table § below) makes it 
paramount to test the approximation schemes of the discussed papers and the hypotheses 
about the importance of various contributions. The results of our analysis are suitable for 
this purpose since the phenomenological amplitude determined by fittings is consistent at 
least with the treated data. 
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The following Table ^ contains the listing of (nonzero) contributions to the threshold 
quantities Di and D2 found in the solution with Xbf = 1-16 (see the variant "DRNbior" of 
the Table |^ for values of scattering lengths). We also quote here the predictions of papers 
T^ , |l^] and the results ALL and TRI of the linear fit from the Table ^ for an easy 



comparison. The underestimate of the D2 value by the linear fits ALL and TRI is due to the 
pointed above tangible departure of the { — \-n} channel's solution from the linear pattern. 

It should be noted once more that there is no absolute meaning of the separate contri- 
butions because of the field redefinition freedom. However, the "on-mass-shell" parameters 
9oy 9ii 92 and (73 of the Air vertex are stable and the tttt scattering lengths are well defined in 
our approach though the "off-shell" contributions of the parameters go, gi, g2 and 5^3 in the 
Table ^ are model dependent. Hence, the most general inferences which might be derived 
from this Table are the following: 

1. The resulting values of the threshold amplitudes and the quantities Di and D2 are 
rather small differences of large contributions from various mechanisms. 

2. In both quantities the OPE mechanism meets the strong competition from all the rest 
ones, the A being of importance for Di while the iV^. — for D2- 

3. Within the overall OPE contribution the influence of the D-wave parameters g2i gs 
reaches 30% in both quantities Di and D2. 



5.3 Conclusion 

We have already seen in the previous subsection that in the experimental side the resolution 
of an ambiguity in two threshold amplitudes and upon which the Olsson-Turner 
approach is relying is the matter of the precision of experimental data on 7r27r total cross 
sections. We also noticed the evidence of the inapplicability of the linear fit for the { — \-n} 
channel — this enlarges the systematic errors of and values. Here, one can see that 
there are also unknown systematic errors in the theoretical field. To fill in this gap the large 
amount of experimental and phenomenological information on isobar physics is necessary. 
Therefore, the approach is losing the advantage of simplicity which had been making it so 
attractive. 
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A2 
A, 
A, 
A, 
A12 

^13 
^15 

Aid 
An 



310.80 
194.66 
531.15 
597.88 
-35.15 
-69.56 
-312.40 
-159.25 
-186.93 
-170.52 



± 
± 
± 
± 
± 
± 
± 
± 
± 
± 



35.44 

6.66 
13.61 
21.26 
95.36 
8.05 
7.48 
4.08 
16.28 
21.80 



310.80 
194.66 
-1062.31 

1195.77 
-35.15 
-69.56 
624.80 
318.50 

-186.93 
341.04 



£2 
± 
± 
± 
± 
± 
± 
± 
± 
± 
± 



35.44 
6.66 
27.23 
42.52 
95.36 
8.05 
14.96 
8.17 
16.28 
43.60 



Ol 

02 
03 
04 



102.14 ± 

488.02 ± 

-101.07 ± 

-161.90 ± 



251.12 
31.38 
33.53 
32.09 



102.14 ± 

-976.04 ± 

-101.07 ± 

323.80 ± 



251.12 
62.76 
33.53 
64.19 



r2 



300.89 ± 



7.75 



-601.78 ± 



15.50 



D2 



-46.75 
-238.69 
93.27 
-14.91 
5.88 



± 
± 
± 
± 
± 



15.48 
51.93 
23.77 
66.48 
11.12 



93.50 
42.32 
-186.54 
-115.75 
-19.53 



± 
± 
± 
± 
± 



30.96 
55.87 
47.53 
133.00 
33.28 



Ri 
R3 

Re 
R7 



52.87 ± 

65.55 ± 

0.46 ± 

2.41 ± 

-7.91 ± 



108.57 
26.58 
12.95 

211.45 
4.73 



-1079.45 ± 

-131.09 ± 

-90.90 ± 

28.60 ± 

-40.14 ± 



192.22 
53.15 
19.91 

505.56 
48.72 



N2 



128.70 ± 
152.34 ± 
-12.29 ± 



32.77 
45.49 
22.62 



2.81 ± 
-304.69 ± 
15.16 ± 



22.75 
90.97 
21.78 



Sum 



313.92 ± 



101.37 



-1407.03 ± 



229.23 



m 
m 
m 



339.69 ± 
344.89 ± 
334.48 ± 



28.63 
31.24 
13.01 



-1186.96 ± 
-1179.15 ± 
-1231.21 ± 



123.64 
136.66 
7.81 



ALL 
TRI 



327.78 ± 
300.46 ± 



14.90 
14.90 



-1025.54 ± 
-1008.16 ± 



29.80 
29.80 



Table 8: Contributions to the values of threshold Di and D2 (in [GeV] ^). 
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6 Chew— Low Extrapolation 



The existence of the pole at r = /x^ in the OPE contribution is the keystone of the Chew- 
Low approach. However, there are several routines for the extrapolation to this point — 
the relevant discussion on the applications of the approach might be found in the review 
paper (see also |58[ for a phenomenological introduction); the tests of some variants 



performed long time ago were reported in the paper [g. The amplitudes provided by our 
fits present the possibility to test the approach by modeling the experimental data and to 
arrive at conclusions which are independent of the finite precision of the input. 

6.1 Extrapolation Function 

From the very beginning it has become evident that the extrapolation function constructed 
in terms of the total cross section could not provide a base neither for the linear extrapolation 
nor for the quadratic one (with the only exception for the case of the 7t2tt amplitude build 
of a constant OPE term plus a constant in the same spinor structure S). Therefore, to 
save the space we discuss only the extrapolation function Fm{j) defined in terms of the 
quasi-amplitude: 

Fm r = — X , — . 36 

The principal feature of this extrapolation function is the one, reflecting the vanishing of the 
cross section at r = 0. It will be commented later on during the discussion of results. 

Since — according to the general idea of the Chew-Low extrapolation — only the OPE 
term of the i^l-n amplitude contributes to the quantity F/vf (/i^) it is convenient to define the 
auxiliary function 

^"^^^"^"^J = (-r)(2,.^^)2 ^ MX) " // dB.dv, ■ 

n(r) 

In the above equations both the LHS and the integrals of the RHS depend on the reaction 
energy s and the dipion invariant mass St^tt which are assumed to be fixed in the course of the 
extrapolation; the amplitude Mqpe is obtained from M by setting to zero all contributions 
but the OPE one. 

The analytic calculation of the function Fm{t) is possible only for the case of the simplest 
models of the amplitude. Therefore, we calculate the function Fm{t) obtained in various 
solutions for our phenomenological amplitude numerically. However, because of the collapse 
of the integration domain f2(r) outside the phase space of the 7r27r reaction the numerical 
calculation of the function Fm{t) at r = /x^ also becomes impossible. (The Monte-Carlo 
based utilities of high energy physics can not generate events outside the phase space for the 
numerical integration.) 
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The integration in the RHS of eq. ( pT] ) results in the rational function of r and s-,^.^. 
Neglecting the imaginary part of the OPE amplitude this function might be cast as the 



quadratic form in the OPE parameters go, gi, g2, gz of the paper 



(go) 




(90 \ 


9i 


(^) 


9i 


92 




92 


V 93 ) 




\ 93 J 



OPEJ 



where the upper triangle of the symmetric matrix $ is explicitly given by 
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(38) 



'^'00 

$11 

$22 
$33 



2; %i = en; $02 = 2(991 + A,)/9 ; $03 = (90| - Ai)/9 ; 
(3e^ + Ai)/6; $12 = i9el + A^)eR/9■, $13 = (9^^^ - 7Ai)^^^/18 
2{A5ej, + mlA, + Al)/A5; $23 = (45^^ - A?)/45 ; 
(454 + 50elA, + Al)/90 ; 

~ [(^i; - - 9g..r/4] ; 9^ = {r + 3fi' - 3s^^) /4 . 



(39) 



The discussed function (|37| ) is well defined outside the reaction phase space and at r = 
gives the value in question -Fope(/^^) = -^m(/^^)- 

In fact, the actual calculations had been being performed with the function -Fope(t) 



-^ope(^) 



[FoPEir)y 



;(40) 



TTTV '—^^iriv 



since the experimental (and the simulated) data are being presented for a strip s^^J — Astttt < 
Stttt < + AsvTTT in the Chew-Low plot in the (r, 3^^^) variables with the non-negligible 
width 2As7r7r. The analytic calculations result in much more complicated expressions than 
that of eqs. (0), so we do not present the final answers here. Because of the rapid growth 
of the TTTT amplitude with Stttt at the threshold the difference of the simple function (|37| ) with 
the above one (|40|) was found to be reaching 20% at As^^^ = 0.15/i^. 

The nonnegligent spread in s,r7r has another important issue for the extrapolation. The 
specific bin (a) is then characterized by the rectangle (r^"-* ± Ar, s^J ± As,r7r) in the Chew- 
Low plane (r, s,r7r)- The cross section for the bin is given by the integral 



^w(l|M|r;r) 



„MAX 

7V7V ^^TTTT 



ds^^ R{\\M\\' 



r) 



(41) 



where i?(||M||^; Stttt, t) stands for the matrix element integrated over the rest 2 variables 
(which include the vrvr scattering angle). The upper limit 



MAX 



= MAX{4"j + As 



(r)} 



(42) 



is independent of r only for bins (a) for which the strip s^^ — Asnn < s^n < s''^ + Asnn is 
going strictly inside the physical domain of the Chew-Low plot. Therefore, in the case of 
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bins located at the boundaries of the physical interval [r (s),r+(s)] the physical space in 
the Stttt variable is bounded not by the value Stttt + Astttt independent of r but by the curve 

s+^(r) = ^{^{s + m^- /i^) + 2m^n^ + ^t{t - Am?){s - (m + - (m - ^if)] . (43) 

As a result, the value of the phase space (T(Q,)(l;r) for such bins does depend on r as 
well as integrals of powers (s,r7r)" do. This makes necessary to withdraw such bins from the 
extrapolation base regardless of the kind of the extrapolation function. For example, both 
the cross section (T(q,)(||M|P; r) and the quasi-amphtude have the breaking points at 2 values 
of r — at the ones for which s5jr^(r) = St^tt + As^^r. This phenomenon is clearly seen in the 
Fig. 10. In the course of a practical data treatment the selection of bins is an easy problem 
which is solved by calculating the empty phase space for the considered array of bins and 
keeping on the ones with the constant value of the phase space. 

Would one know in advance the s,r7r dependence of the partially integrated matrix element 
R{^s-n-ni t) it would be possible to make corrections for the bins intersecting the boundary curve 
s+^(r) and to include more points into the Chew-Low extrapolation. Our curves in Fig. 10 
had been corrected by the empty phase space — evidently this is insufficient. (Since our 
solutions fit well the s-^t^ distributions which are very distinct from the phase space there are 
no much room for wondering.) Another possibility to enlarge the base of the extrapolation 
by cutting more narrow strips in s,r7r depends completely on the experimental statistics. 



6.2 Simulations of Chew— Low Extrapolation 

The statistics of our data (see experimental points in Fig. 10) can not provide a confidence for 
the results of a real-data extrapolation. Therefore, we were simulating the distributions with 
the help of the theoretical amplitudes and constructing the extrapolation function Fm{t) 
for a reasonable number of bins. The input data errors then become negligible, hence, the 
problem might be investigated in its pure state. 

Here, we discuss the simulations of r distributions (for the fixed strip in 5^,^) performed for 
several energies, namely, for PLab = 335, 420 and 460 MeV/c. The extrapolation functions 
calculated for three types of theoretical amplitudes at PLab = 460 MeV/c are shown in 
Fig. 10. The shown data are simulated with the binning and the precision which are only 
computer-dependent; the same binning of the available experimental data suffers from a lack 
of statistics — this is clearly demonstrated by the empty experimental bins in the discussed 
pictures. 

The simulated data were subject to extrapolation to the point r = /i^. The true limiting 
value for each amplitude is being calculated with the use of eq. (^01). The linear {lin) and 
the quadratic {squ) extrapolation patterns are selected for demonstrations. The results of 
extrapolations are collected in the Table ^ where the fulcrum numbers of the true limiting 
values specific for the considered amplitudes are given in the bottom boxes. We display here 
the variation of the extrapolated values with the choice of the left bound xi, the right bound 
T2 being kept fixed. 

The o columns of the Table ^ are the undoubted grounds for the crucial inference that 
even in the simplified case of the pure OPE mechanism the linear extrapolation method 
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PLab = 335 MeV/c, = 4.15//2 



n 






Ohn 


Osqu 


giin 


9squ 






17 


-5.519 


-1.155 


0.6111 


0.7889 


1.201 


0.866 


0.1612 


0.3578 


16 


-5.262 


-1.155 


0.6204 


0.7883 


1.176 


0.909 


0.1712 


0.3582 


14 


-4.749 


-1.155 


0.6382 


0.7871 


1.138 


0.987 


0.1898 


0.3669 


Iz 


-4.235 


-1.155 


0.6546 


0.7892 


1.112 


1.067 


0.2087 


0.3723 


10 


-3.722 


-1.155 


0.6709 


0.7875 


1.100 


1.140 


0.2286 


0.3697 


8 


-3.208 


-1.155 


0.6858 


0.7864 


1.099 


1.223 


0.2479 


0.3589 


6 


-2.695 


-1.155 


0.6998 


0.7906 


1.111 


1.341 


0.2607 


0.3953 










0.7870 






0.7870 






0.6036 






0.6036 






0.0000 






0.0000 





PLab = 420 MeV/c, s^^ = 4.15/x2 



n 


ri//i2 




Olin 


Osqu 


Qlin 


9squ 






17 


-10.01 


-1.283 


0.4497 


0.7097 


0.4011 


0.3829 


0.4971 


0.5755 


16 


-9.497 


-1.283 


0.4672 


0.7090 


0.3848 


0.4435 


0.4971 


0.5965 


14 


-8.470 


-1.283 


0.4999 


0.7098 


0.3715 


0.5568 


0.5051 


0.6259 


12 


-7.444 


-1.283 


0.5315 


0.7055 


0.3838 


0.6502 


0.5208 


0.6395 


10 


-6.417 


-1.283 


0.5598 


0.7028 


0.4154 


0.7361 


0.5384 


0.6496 


8 


-5.390 


-1.283 


0.5845 


0.7070 


0.4632 


0.8261 


0.5561 


0.6688 


6 


-4.364 


-1.283 


0.6089 


0.6990 


0.5321 


0.8703 


0.5782 


0.6715 










0.7054 






0.7054 






0.4751 






0.4751 






0.0000 






0.0000 





PLab = 460 MeV/c, s^^ = 4.45//^ 



n 








^squ 


9lin 


9squ 




•^squ 


21 


-12.58 


-1.797 


0.1849 


0.4911 


-0.0071 


0.908 


0.2755 


1.485 


20 


-12.06 


-1.797 


0.2015 


0.4916 


0.0174 


1.006 


0.3477 


1.462 


18 


-11.04 


-1.797 


0.2337 


0.4917 


0.8569 


1.200 


0.4837 


1.405 


16 


-10.01 


-1.797 


0.2643 


0.4917 


0.1830 


1.378 


0.6084 


1.325 


14 


-8.984 


-1.797 


0.2928 


0.4933 


0.3074 


1.535 


0.7152 


1.230 


12 


-7.957 


-1.797 


0.3198 


0.4951 


0.4568 


1.658 


0.7997 


1.131 


10 


-6.930 


-1.797 


0.3459 


0.4917 


0.6242 


1.737 


0.8584 


1.042 


8 


-5.904 


-1.797 


0.3695 


0.4895 


0.8043 


1.726 


0.8938 


0.969 


6 


-4.877 


-1.797 


0.3906 


0.4939 


0.9711 


1.653 


0.9096 


0.922 










0.4886 






0.4886 






0.3591 






0.3591 






0.0000 






0.0000 





Table 9: Results lin (squ) of the linear (quadratic) Chew-Low extrapolation to the point 
r = for the varying left bound ti. The theoretical amplitudes used for data simulations 
correspond to: o — the solution with OPE contribution only; g — the solution with all 
mechanisms; x — the solution with all mechanisms excluding OPE. The numbers given in 
the bottom boxes show the true values at r = //^. 
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generally underestimates the value in question and results in the inappropriate systematic 
error of 25-35%. The impression of some improvement with the shift to the extreme right 
position of the extrapolation database in the r interval is misleading since in the conditions 
of the reduced extrapolation base the effect of nonzero errors of the real experimental data 
must make the result even more ambiguous. We must also note that the OPE amplitude 
(i.e., 6) does not fit at all the overall data (see Table ^ for values of x^)- Besides, it was 
several times pointed out that only the "gauge-covariant" set of contributions makes sense 
due to the field-redefinition freedom. 

In the more realistic case g when all mechanisms are being present (the amplitude g of 
the solution with Xtif ~ ^-^^ is at least compatible with the overall database) there are 
no advantages in both the linear and the quadratic extrapolation methods; the coincidence 
of the results with the exact answers seems to be of a random nature. The 200-300% 
deviation makes it unreliable to use both the linear and the quadratic extrapolations even 
for estimations. 

What is really disappointing it is the examination of the columns x. The nontrivial 
answers in this case raise suspicions that the extrapolations follow the dictate of the experi- 
mental data rather than the theoretical amplitude. Indeed, the theoretical amplitude x fits 
well the data but it has no pole at t = fi"^ at all! 



6.3 Discussion 

In view of the negative general conclusion on the applicability of the Chew-Low extrapolation 
approach at the considered energies we need to shed more light on its origin. 

The separate clarification is necessary for the case of the pure OPE mechanism since the 
application of the Chew-Low approach is based on the hypothesis of the OPE dominance. 
Indeed, if the nature follows the simplest pattern of the OPE dominance in the 7r27r reaction 
then 6-8 data points are enough for successful extrapolation — the small discrepancy dis- 
played in the Table |^ (which must be growing with the energy) is due to different accounting 
of isospin breaking in the main program and in the extrapolation function (0). 

It is not so difficult to realize that the small departure of the extrapolation function ([37| ) 
from the linear pattern is solely due to the participation of the D-wave parameters g2 and 
gs — see the quantities $22; '^'23, $33 given by eqs. ([391). In the absence of the latter the 



internal integrations of the leading order ChPT amplitude in the formula (^) result in the 
linear function of r. At the considered energies the influence of the quoted parameters on 
the TTTT amplitude itself is negligible. The deviation of the extrapolation function from the 
linear shape is small but its effect on the results of extrapolations is found to be drastic. 

The off-shell appearance of the tttc amplitude in the 7r27r reaction acts the part of the 
magnification lens in respect to Z)-wave parameters — we have already seen this in the pre- 
vious section when discussing their contributions to the threshold amplitudes. However, the 
above phenomenon does not present an obstacle by itself since the quadratic extrapolation 
for a pure OPE amplitude is proved to be exact and stable. 

It is the complicated form of the 7r27r amplitude revealed by our data fittings which rules 
out the possibility of a reliable application of the Chew-Low extrapolation in the simplest 
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manner. This conclusion is derived in terms of a particular ansatz of the extrapolation 



function ([36|). Let us now discuss why the conclusion is of more general nature. 

There is the difference of our function Fm{t) with the ones defined in terms of cross 
sections — we are extracting the square root of the eq. (|36|). The results of the extrapolation 



of the square of the function Fm{t) in terms of the ansatz 



l2 



[Fm{t)Y = [Fo + F^{T-^') + F2{T~^y\ (44) 

are completely equivalent to the ones displayed in the Table ^ At the same time the 
extrapolation via 

[FM{T)f = h + h{r~^?) + F2{T-^?f (45) 

provides worse results. 

Hence, we see no reason in keeping on the cross-section form. It is not the point which is 
capable to disapprove our conclusions. So let us discuss another feature which is implemented 
into our extrapolation function (|5BD. 

There is the property of the pure OPE cross section o"|r^o = which was displayed long 
time ago as by the nonrelativistic calculations (see, for example the textbook by Kallen 
p9| ) as well as by the relativistic ones (like that of the paper by Naisse and Reignier). 
The work P| by Baton, Laurens and Reignier provided the phenomenological test of this 
property by the high energy (pLab = 2.77 GeV/c) data — since then it was being built into 
the applications of the Chew-Low procedure as a standard feature. Nevertheless, the known 
failures of applications of the Chew-Low approach were associated with the relying on the 



very property we are discussing here — see the review |Q by Leksin. The more close analysis 
shows that in some cases the foothold on the property (J\t^o = in the definition (|SB| ) is the 
reason of overshooting of the quadratic extrapolations which are accurately following the 
extrapolated data in the physical region. 

Let us now briefly remind what is the theoretical status of the hypothesis that air^o = 0. 
Definitely, it is the exact property of the pure OPE mechanism. What the value of the matrix 
element (7r27r3A^(g)|S'|7riA^(p)) at p = g is in the general case is the kinematical problem in 
part. (It should be noted that the point for which p = q is located outside the physical 
region both for the t^N nirN reaction and for the 47r vertex since the condition p = q 
implies Stttt = A*^-) 

The 5* structure of the amplitude (0) gets the same multiplier (— r) in the unpolarized 
matrix element (|l^) as the OPE contribution does. This structure gives rise to spin-flip am- 
plitudes which are the only amplitudes of the considered reaction surviving at the threshold. 
However, besides this spinor structure there are three more; kinematically, their contribu- 
tions to the quantity ([l3| ) are determined by the matrix ([T^) . 

Thus, to make the matrix element (|1^) vanish at p = q all entries of the matrix ( p^ must 



become zero simultaneously. The analysis of the explicit expressions (for which we have no 
room here) shows that three conditions are necessary: 1) s = {m + a/s^)^; 2) coUinear final 
pions k2 = ks] 3) Chiral limit fi = 0. The last condition is also the only general reason to 
make the considered structures vanish dynamically. 
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Therefore, in the real dynamics the quantity ||M||^_^q {a{T = 0)) stands for expressing 
the Chiral symmetry breaking which is similar to the vrA^-elastic S term. It depends on the 
energy but seems to be rather small. Nevertheless, it prevents to make the safe simplifica- 
tion in the definition of the extrapolation function, namely, to divide the quasi-amplitude 
by a/— r in our case. Thus, we arrive at the conclusion that, from one side, the complicated 
dependence of the physical amplitude on r makes useless the linear and the quadratic extrap- 
olation methods even in terms of the quasi-amplitude, from the other side — the presence 
of the Chiral symmetry breaking in the true amplitude forbids to soften the dependence in 
the extrapolation ansatz. 



6.4 Remark on D— Wave Parameters 

The parameters g2, of the cross-symmetric ansatz of the tttt amplitude (see [^], [21 



determine the values of the D-wave scattering lengths, therefore, we call them the 
Z^-wave parameters here. Because of the crossing symmetry the same D-wave parameters 



determine also the slopes of the S'-wave amplitudes (see eqs. (59), (60) of the ref. and 
eqs. (73)-(104) of the ref. |23|). 

Due to the considerable growth of the 1 = amplitude the contribution of the slopes to 
the integrated over St^tt cross section (extrapolated by the Chew-Low method) appears to 
be large. 

This very phenomenon makes it possible to extract the discussed parameters from the 
near-threshold 7r27r experiment, otherwise it is simple to verify that at the considered energies 
the contribution to the hit and 7r27r cross sections by the D waves themselves is negligible. 
This also presents the additional motivation for the careful handling of the crossing properties 
of the amplitudes of tttt and 7iN —>■ tttcN reactions. 



7 Conclusions 

Throughout the paper we were making the inferable statements along the discussions. Here, 
we remind the most important ones and develop the general conclusions. 

7.1 Data 

Our present work is devoted to the analysis of the near-threshold data on the tt2tt reaction. 
The data base described in sect. 3. consists of the experimental total cross sections and 1- 
dimensional distributions. The full-kinematics data of the work also had been presented 
in the same form. This needs some comments. 

The available 1023 full-kinematics events of the quoted work constitute the solid ground 
for the total cross section; 10-14 bins of a 1-dimensional distribution have a good filling with 
the averaged number of 60-100 events per bin; the filling of 8 x 8 bins of a 2-dimensional 
distribution is satisfactory (15 events per bin in the average) while the filling of 6 x 6 x 6 of the 
3-dimensional ones and 4x4x4x4 of the 4-dimensional bins is poor. In this conditions of 
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the difficult choice between the poor filhng of bins and the loss of the kinematical information 
we formed multiple 1-dimensional projections for the data to bring to light the behavior in 
the crucial variables. 

We created some additional ID projections of the data fl^ and tested them with the 
obtained solutions — the description was found to be excellent. 

The use of the numerous lower dimensional distributions acts the part of a kind of the 
tomography method. The fittings showed that at a distance from resonance poles this works. 
The improved statistics of the contemporary experiment |7^, ^ definitely must make 

the direct use of the full-kinematics data more preferable. 

The general properties of the treated data base are found to be: 

1. The precision of data is insufficient to improve the accuracy of the determination of 
characteristics of the vrvr scattering. However, this is the problem not only of the data. 

2. The coverage of the behavior of the low energy amplitude is good; practically, there 
were no losses of convergence in the course of fittings. The numerous solutions reflect the 
complexity of the reaction amplitude — the current experimental setup can not be charged 
for this (we shall continue this discussion below in the subsects. 7.2., 7.3.). 

3. The distribution data of the {+ + n} channel are found to be extremely important 
for the resolution of parameter correlations. In the absence of these data the convergence 
of fittings becomes tremendously slow, the number of iterations being increased by several 
orders. 

To resume we state that the considered data base is in principle sufficient for determi- 
nation of the phenomenological near-threshold 7r27r amplitude. Certainly, the need in the 
better quality of data in respect to the precision of the obtained parameters and the mul- 
tiplicity of solutions is obvious. First, the quality of the newest data will be much better 
and, second, it seems oversimplified to charge only the data with this problem. The more 
detailed discussion of its origin will be given in the next subsection. 



7.2 Amplitude and Major Results 

Our amplitude is built on the rather conservative basis. However, the orientation towards 
the modern ChPT approach in constructing the model might be lacking of an intrinsic tool 
for making the right explanation, if failed, whether an inconsistency of data or the neglect 
of higher order terms are responsible for the inappropriate fit. If successful, the approach 
provides the only conclusion that at the considered order ChPT is compatible with data — 
a substantial estimate of the systematic error of determination of the low energy constants 
is impossible. 

The approach of HBChPT deserves the separate remark. Unlike the case of the vrA^- 
elastic scattering this approach is, probably, inapplicable to the case of the vrA^ ttttN 
reaction considered here. One should remind the basic keystones of HBChPT: 1. Nonrela- 
tivistic limit; 2. Small-pion-momentum expansion; 3. Heavy baryon approximation. Then, 
let us consider the identity 

u{q){ki -h- = -2mu{q)i'j5u{p) , (46) 
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which is specific to the relativistic form of the 7r27r amphtude. The identity can not support 
the above points 2., 3. simultaneously. 

In the present paper we consider the amphtude of the 7r27r reaction built of numerous 
resonance contributions (including the separately treated OPE mechanism) and the smooth 
polynomial background (see sect. 2.). Its complicated appearance reflects the influence of 
the (generally, off-shell) processes like tttt tttt, ttN ttN, ttN tcN^, tcN ttA on the 
near-threshold region of the discussed reaction. 

The fittings strongly confirmed the importance of all quoted exchange mechanisms. This 
result can not be considered as the totally new one. For example, the importance of isobars 
already had been stressed in the paper in terms of the sophisticated analysis of the 
amplitude form. 

The near-threshold region 280 < PLab ^ 500 MeV/c can not be considered as the 
selfcontained one because of large A, A^* widths. The isobars extend their strong influence 
up to the very n2n threshold (this is demonstrated by the Table In the absence of isobar 
contributions the considerable improvement of the fit due to the imaginary background serves 
as an indirect evidence of the importance of the discussed isobar mechanisms (see Table 
The interrelation of these mechanisms with the OPE one will be discussed below. 

In view of the discussion of the role of background parameters (see subsect. 2.4.) it is 
not so surprising that only few of them are found important in the fittings when all exchange 
mechanisms are being present. The fact that the parameters of the isospin amplitude D are 
found to be consistent with zero reflects the negligible influence of higher r-resonances (like 
SE/vr7v(^^)) on the amplitude in the considered energy region. Even the nonzero contributions 
to this isospin amplitude from isobar exchanges almost cancel each other — this might be 
derived from an approximate equality of the resulting theoretical cross sections of {—Op} 
and {+0p} channels. 

The OPE contribution is in the center of our investigations. In all variants the improve- 
ment of Xdf with the inclusion of OPE is found to be statistically important (see Tables ^ 

The 47r vertex of the OPE graph is taken in the direct amplitude form which contains 
4 parameters go, gi, g2 and gs and respects the isospin, crossing and approximate-unitarity 
properties of the tttt amplitude off the mass shell. In terms of the vrvr-scattering lengths the 
values of these parameters in the best physical solutions are: 



Xdf 


1=0 

"o 


"o 








1.161 
1.203 
1.205 
1.212 


0.07±0.12 
0.07±0.11 
0.17±0.12 
0.19±0.11 


-0.056±0.036 
-0.076±0.073 
-0.053±0.039 
-0.059±0.045 


0.045±0.017 
0.047±0.025 
0.054±0.022 
0.053±0.022 


0.0052±0.0031 
0.0023±0.0027 
0.0054±0.0025 
0.0062±0.0027 


-0.0005±0.0014 
-0.0013±0.0014 
-0.0002±0.0012 
0.0006±0.0015 


m 


0.26±0.05 


-0.028±0.012 


0.038±0.002 


0.0017±0.0003 


0.0001±0.0003 



Here, we list also the currently adopted experimental values of the compilation pO[| . 

Generally, the solutions display that the precision of determination of the Z)-wave pa- 
rameters g2 and g^ is not worse than that of 5^0, gi parameters. This was attributed to the 
characteristic energy dependence of the isospin-zero S wave (see subsect. 6.4.). The param- 
eters of the latter via crossing and kinematics are connected to Z^-wave scattering lengths 
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— this is clearly demonstrated by the Oj^" errors in the above list. 

The poor precision of determination of the vrvr-interaction parameters has its origin in 3 
principal reasons: 1) the data accuracy; 2) the competition of other (nonOPE) mechanisms; 
3) the incomplete nature of the current experimental setup. 

The first point is evident — the data accuracy will be improved soon (look |^ for the 
survey of experiments). 

The second one stems mainly from the isobar exchanges — this is clearly revealed by 
fittings: whenever an isobar mechanism is absent it is easy to find a solution with ag^° in the 
range 0.20 — 0.30; it drops to the 0.00 — 0.20 range if all isobar mechanisms are involved 
(it is interesting that the effect of p exchanges is quite opposite). In other words, the larger 
values of vrvr-scattering lengths are gained when OPE is forced to stay for an essential but 
missed isobar contribution. 

The important question is about the nature of parameters which correlations with the 
OPE set go, gi, g2, gs are so devastating. Unfortunately, these are not the parameters of 
the DE-type graphs which might be estimated from the decay characteristics or might be 
known (like the g-j^NN constant) from the low-energy ttN phenomenology. Instead, the root 
correlations are due to parameters Ri-R^ and D1-D4 of the SE^Ar(A^*), SE^Ar(A) graphs (see 
Tables ^, 0, The explanation is simple. Parameters of the graphs of the SE type play the 
same role in respect to the ones of the DE type as the background parameters Ai-Aig (and 
their analogs iig-ise) do in respect to all exchange parameters. Indeed, outside the resonance 
region the contraction of any pole in the DE graph leads to the single-pole contribution, i.e. 
to the one described by the SE graphs. 

Here, one observes the interrelation of the point 2. with the first one since the data 
restricted to the region below resonances can not help much to fix up parameters of processes 
like 7iN ttN^, rnrN, nN ttA, etc. 

We found the numerous set of solutions describing the data at the acceptable level of x^- 
The origin of this phenomenon must be explained in part by the importance of all 4 spin 
structures of the 7r27r amplitude (see eqs. (^), (|T0|)). This is the reason number 3 for the 
poor accuracy of our final results. 

Indeed, the unpolarized data measure only one combination of spin structures (namely, 
the combination given by the matrix element (|l5)) leaving them almost free to stand one 
for another. Therefore, the absence of polarized measurements in the energy region PLab ~ 
500 MeV/ c and the abundance of mechanisms specific to the considered reaction is the reason 
of huge correlations of OPE parameters with the rest ones on the available data base. 

Meanwhile, the extreme importance of the nucleon spin in the considered reaction at 
higher energies had been recently reported by Svec in the paper [jTSf. In the results of 



our analysis and modeling the Chew-Low extrapolation function we see the clear signal of 
nontrivial spinor structures. This claims for the polarization measurements of ttN tcttN 
reactions at the discussed energies. Up to now the known polarization measurements of 
the 7r27r reactions had been performed at considerably higher energies, for example, at 5.98 
GeV/c and 11.85 GeV/c |51[ and at 17.2 GeV/c |3|. Their analyses 0, g already proved 



such measurements to be detailed sources of information on the tttt interaction (at high 
energies) . 
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Certainly, the complete polarization experiment requires the analysis of the polarization 
of the final nucleon — in the near future this is hardly to be carried out for such rare 
processes as the considered one. Nevertheless, the examination of the spinor structure of the 
considered amplitude (^) displays that the almost exhaustive information might be obtained 
already from the experiment with the polarized target. Indeed, there are two independent 
asymetries 

a(s; ki) + (t(s; -ki) ' 

where s is the vector of the nucleon polarization and k^ are the projections of vectors 
k-i- = k2 ± ka to the plane which is orthogonal to s. Their measurements must provide an 
information on two additional combinations of four spinor structures of the decomposition 
(^) which are independent from the combination of the matrix element (|T3p. 

7.3 General Conclusion 



In view of the results of modeling the Chew-Low extrapolation and the Olsson-Turner ap- 
proach we conclude that only the approach based on the extensive phenomenological model 
(a la Vicente- Vacas) can help in investigations of the considered reaction. The investigations 
require to develop the common analyses of related processes like nN — >■ irN^, nN — > vrA — 
i.e. the processes described by the Lagrangian terms listed in Tables ||, | — and, hence, to 
extend the energy region up to pLab ~ 1 GeV/c. In other words, the problem of determi- 
nation of the low-energy Trvr-scattering characteristics is a part of the more comprehensive 
problem of investigation of tt2tt dynamics at low and intermediate energies. 



The results of contemporary experiments [40, 73, 74| must provide much more precise 



determination of vrvr parameters since their correlations with the unknown parameters of 
isobar mechanisms originating from the discussed above spin structures will be resolved in 
part due to the improved data accuracy. In this respect the role of the polarized data is 
difficult to over-estimate. 

The amplitude which structure will be fixed by the analysis of statistically significant 
data might gain the wide range of applications beyond the testing of ChPT predictions. For 
example, it might be used for the common analysis of vrA^ — >■ ttttN, imN and vriV- 

elastic data, for investigations of the rj production in the process vrA^ rjN, for correcting 
experimental distributions obtained at devices with the restricted geometry and for other 
investigations at intermediate energies. This is ensured by the fundamental role of the 
irN ttttN reaction in nuclear and particle physics. 
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Figure 2: Experimental points of total cross sections and the theoretical 
curve for the best physical solution (quasi-amplitude (M) in GeV~^). 
Triangle points were excluded from fittings. 
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Figure 3: Experimental distributions for the {+ + n} channel from the 
paper by Kirz and theoretical curves (normalized quasi-amplitude 

(M)norm). 




Figure 4: Experimental distributions for the { \- n} and { \- n} 

channels from the paper by Jones and theoretical curves (normalized 
quasi-amplitude (M)norm)- 




Figure 5: Experimental distributions for the { — \- n} channel from the 
paper by Kirz and theoretical curves (normalized quasi-amplitude 

(M)norm). 
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Figure 6: Some of experimental distributions for the { — \- n} channel 
from the paper by Blokhintseva and theoretical curves (normalized 
quasi-amplitude (M)norm)- 
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Figure 7: Experimental distributions and theoretical curves for angular 
spectra of the { — \-n} channel from the Saxon's paper 0] and the same 
spectra build of the Blokhintseva data [|T^ (normalized quasi-amplitude 
(M)norm)- The Saxou data had not been used in fittings. 




Figure 8: Experimental distributions and theoretical predictions for 
spectra of the { — \- n} channel from the paper |j7l| by Saxon and from 
the paper by Kirz (normalized quasi-amplitude (M)norm)- 



2000 



> 



A 
V 




2000 



2000 



0) 



V 



1000 



500 





TT** tt" n 






" 1 1 


1 1 1 1 1 1 1 



0.05 0,1 o,i; 

T (Gev) 



V 



2000 
1500 
1000 
500 












ra 




1 1 1 1 1 1 1 1 1 1 



0.05 0.1 0.15 

T (Gev) 



V 




2000 



Q 1 500 



1000 



V 




500 — 



800 
P, 600 - 
400 — 
200 — 




> 



V 



square tt tt | 
I dot re* n° p T il 




Figure 9: Linear fit of total cross sections of all five channels with 7 
parameters (quasi-amplitude (M) in GeV~^). 
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Figure 10: Simulations of extrapolation data (Fm)'- dot — the amplitude 
of the best physical solution; square — the amplitude without the OPE 
contribution; circle — the pure OPE amplitude; cross — the available 
experimental data of the paper 
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